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DIURNAL VARIATIONS IN THE NUMBER OF SHOWER METEORS 
DETECTED BY THE FORWARD-SCATTERING OF RADIO WAVES 


PART I. THEORY! 


By C. O. HINEs 


ABSTRACT 


Various observational factors affect the number of shower meteors which can 
be detected, in a given time interval, by the scattering of radio waves from the 
ionized meteor trails. The pertinent factors in the case of forward-scattering, 
where the transmitter and receiver are widely separated, are expressed approxi- 
mately as functions of the position of the shower radiant. In combination, they 
provide an observational weight factor which may change appreciably as the 
radiant moves in the course of a day. The consequent diurnal variation in the 
occurrence of scattered signals may then be determined, and distinguished from 
variations due to random changes in the incidence rate of the meteors. 


INTRODUCTION 


A number of meteor swarms are encountered yearly by the earth, as it follows 
its course round the sun, and these give rise to a succession of meteor showers. 
The number of shower meteors incident on the earth may vary systematically 
as the earth passes through a swarm, and random fluctuations in the incidence 
rate will inevitably occur. The number incident above an observer's horizon 
will show a further variation, in this case diurnal, dependent upon the varying 
elevation of the shower radiant. 

When the ionized meteor trails are detected by the (forward) scattering of 
radio waves transmitted from a distant station, their numbers are often 
found to vary appreciably throughout each day in a manner which cannot 
reasonably be attributed to any of the causes just mentioned. It seems likely, 
then, that the method of observation is itself responsible for a large part of 
the variation observed. 

This possibility is investigated in the present paper, using a theory of 
forward-scattering developed by Eshleman (1). The treatment is approximate, 
and not justified if the transmission path is much shorter than 1000 km.; it is 
illustrated by contour charts for this path length, but is otherwise kept as 
general as possible. The charts are most simply applied only after a trans- 
formation of coordinates determined by the particular position and orientation 
of the path in question, but this stage of the work is given in outline only. 


1Manuscript received March 18, 1956. 
Contribution from the Radio Physics Laboratory of the Defence Research Board, Ottawa, 


Canada. The work was performed under project PCC No. D48-95-11-01. 
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The transformed chart appropriate to one path now in use, and a comparison 
with the corresponding observational data, will be included in a subsequent 
joint communication. 


PRIMARY CONSIDERATIONS 


The numeric count of shower meteors is governed basically by the number 
density in the meteor stream. Random fluctuations in this density may be 
severe, and may override much of the predictable variation introduced by 
observational factors. In calculating these factors then, it seemed reasonable 
to the writer to accept any approximation—even if only qualitatively 
justified—which would give a preliminary indication of the net observational 
effect. With this approach, it has been found useful to calculate many quantities 
at the mid-plane bisecting the transmitter—receiver axis, 7R, and to apply the 
results as representative for other points. This approximation is justified when 
the 7R axis is long enough, since the meteor horizons about T and R then 
overlap only in a narrow band near the mid-plane. (See Fig. 1.) There is, in 
fact, no overlap when 7R exceeds about 2300 km.; but even for TR as short 
as 1600 km., the overlap extends only a third of the way from the mid-point 
M to T and to R. It is thought that the theory will be applicable even when 
the TR axis is as short as 1000 km., since, in practice, the antenna patterns 
will not cover the whole of the overlapping region. The effect of a typical 
antenna pattern will be illustrated in a later section. 
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Fic. 1. Vertical section through the transmitter-receiver axis. 


Forward-scattering from so-called underdense trails may be treated as a 
problem in Fresnel diffraction theory (1), and it is then evident that the 
principal Fresnel zone is potentially the most effective scattering region. 
Taken to an extreme, this fact suggests the following approximation: meteor 
trails which contain the center point P of the associated Fresnel system within 
their lengths are potentially observable, and those that do not are not. The points P 
are properly found as the points at which the meteor tracks become tangent to 
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one or other of the family of prolate spheroids about TJ and R as foci. These 
spheroids approximate to circular cylinders about TR as axis, near the mid- 
plane, and the condition for tangency to cylinders is much more easily obtained. 
If the shower radiant lies in a plane inclined at angle y to the vertical plane 
through 7R, then the points at which the tracks become tangent to the 
cylinders lie in another plane through 7R, this one inclined at angle Y above 
the horizontal direction. (See Fig. 2, which represents projections on a plane 
normal to the plane depicted in Fig. 1, right-bisecting 7R.) The dual role of y 
should be explicitly noted, for it may often lead to confusion: y measures the 
radiant position from a vertical plane, and the tangent points from a horizontal 
plane. 
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Fic. 2. Projection on the mid-plane bisecting the transmitter—receiver axis. 


If the meteor trails are taken to extend between two levels, at heights 
H and # +6 above a spherical earth of radius A(> “> 4), then the 
observationally relevant tangent points will comprise only a part of the 
inclined plane—that part which lies between heights “and “ +46, and which 
lies within the horizon circles about T and R. The slant breadth of this zone is 
given approximately by 6/sin(¥+o), where o is the angle subtended at the 
earth’s center by the line MS, if S is the particular P point which lies in the 
mid-plane at the “ level. The upper and lower sides of the zone approximate 
to straight lines near the mid-plane, parallel to the TR axis. The lower side 
lies at some height / above the horizontal plane through 7R, given by 


{1] h?/sin*y + 2h Bcos Q = #sin?Q4+2A#H +H ?, 


where Q is half the angle subtended at the earth’s center by 7R; if the length 
of TR is 2D, then 
[2] sinQ = D/Z. 


In practice, writing 


[3] H+& —-RoosQ = H, 
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[1] can be approximated successfully by 


[4] h = H(1 — H/2 &cos QP tan) 
for y = 20°, and by 
[5] h = 2H/{1+(1 + 2H/ & cos Q sin’y)}] 


for smaller y’s, without significant error (3%). 

The lower side of the observational zone is separated from the TR axis 
by a distance h/sin y. It extends a distance d on either side of the mid-plane, 
determined by the horizon cutoff and hence by the relation 


[6] d = (F2—h?/sin*y)?—D 
where 
7] P= ((R4H)-AM 


is the radius of a horizon circle at height # . Equation [6] is illustrated by the 
plan view in Fig. 3, although the projections ’ and h cot w appear there in 
place of Aand h/sin y. 


‘| 
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Fic. 3. Projection on the horizontal plane through the transmitter-receiver axis. 


The area of the observational zone can now be determined approximately. 
Its projection on a plane normal to the meteor tracks is reduced by a factor 
sin 8, where 8 is the angle between the plane TSR and the radiant position. 
(8 can be measured in a plane through 7, inclined at angle y co the vertical, 
from the 7R axis as the 8 = 0 direction; together with y, it locates the radiant 
position. See Fig. 8.) The projected area is then given by 


8] A = 2sin B [(Z?—h?/sin*y)*—D] 6/sin(y+c). 


The number of meteors whose trails are potentially observable, in the sense 
stated previously, will be proportional to this area, and so will have the func- 
tional form of the expression [8]. 
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As introduced, 6 is a constant. However, a meteor trail actually extends 
over an undefined distance with varying ionization density, and a weighted 
integration over 6 should therefore be performed in principle. The form of the 
height variation of ionization is, however, effectively independent of the trail 
orientation (cf. Herlofson’s development (3), or Kaiser’s more recent summary 
(4), of the theory of trail formation). The integration would therefore produce 
no new functional dependence on 8, y; it, and the factor 6, may be neglected 
in the subsequent development. 

The height at which maximum ionization occurs does depend on x, the 
zenith angle of the trail axis, but not very critically except when x — 90°. 
The corresponding variation in the appropriate height “ may therefore be 
ignored, and h may be determined from [1]-[5] with ” constant. 


SECONDARY CONSIDERATIONS 


Not only must a meteor pass through the zone of observation for its trail 
to be observed and counted, but also the power scattered to the receiver 
must be sufficient to permit detection. The following functional dependence 
for the power scattered to R has been determined (1, 2): 


[9] Pr « @GrGp sin’a/rir2(r1 +72) (1—cos?8 sin*¢) 


where g is the line density of ionization (assumed constant, but in practice 
at P), Gr and Gp» are the transmitting and receiving antenna gains in the 
directions 7P and RP respectively, r; and rz are the corresponding ranges, 
a is the angle between PR and the electric vector of the (plane polarized) 
wave incident on the trail at P, 2¢ is the ‘forward-scattering angle’ 7PR, 
and 8B is the angle between the trail and the plane TPR. This definition of 8 
is the same as that given previously, when P lies in the mid-plane; it is ap- 
proximately equivalent even at some distance from the mid-plane—so long as 
the cylindrical approximation is valid—and the distinction will be ignored. 

If the receiving antenna accepts only a particular polarization component 
of the wave incident on it, then the sin’a factor in [9] should be replaced by 
cos*y to obtain the functional form of the power actually observed: 


[10] Po « Q’GrGp cos*y/rir2(r71+72) (1 —cos?8 sin’), 


where y is the angle between the electric vector of the wave incident on the 
trail and that component which is accepted by the receiving antenna. The 
modified formula will be adopted in this paper, on the assumption that the 
transmitting antenna propagates, and the receiving antenna accepts, only a 
horizontally polarized component. 

The geometric factors in [10] are easily calculated for the point S. The 
distances 7; and 72 are identical then, and equal to 


(11] r = (D*+h?/sin)}, 
while 
[12] sing = D/r. 


In a Cartesian system corresponding to the polar system defining the shower 
radiant (8, ¥), the zenith at T has direction cosines cos Q, 0, sin Q, and the 
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line TS has direction cosines proportional to h, +hcoty, —D. (See Fig. 4; 
the + alternative accounts for an observation on either side of the TR axis.) 


s 
METEOR TRACK ~ r 





TO CENTRE OF EARTH 


Fic. 4. Geometry of forward-scattering from a trail through S. 


If the electric vector of the wave incident on the trail at S has direction 
cosines X, yu, v, then the condition for a horizontal polarization implies 


[13] AcosQ+yvsinQ =0 
and the condition for a transverse wave implies 
(14] Ah + ph cot y —vD = 0. 


Similar equations apply to the set ’, »’, v’ for the electric vector accepted by 
the receiver; Q and D need only be changed in sign. These relations, together 
with the implicit conditions \?+y?+v? = 1 = )’?+y’2+y’2, can be combined 
to give 
- - 5 (D+h tan Q)’—h’cot’y(1—tan’Q) 
= me ee ee a, a a 
[15] Joos -y] = [AN ban 'ter'| = | ran 0)’ +h cot?y(1+tan’Q) 


~ |D’—h’cot’y|/(D’?+h’cot’y), 
where the approximation is justified over the relevant range 1000 km. S 2D < 
2300 km. 

According to present theory of meteor trail formation (3, 4), the line density 
of ionization at each point along the trail is proportional to the initial mass m 
of the meteor and to cos x. In the mid-plane, this second factor is equal to 
sin 8 cos(¥+c). A further variation with zenith angle results from the varying 
height of maximum ionization, but such effects are small and will again be 
ignored. 

The foregoing results may be combined to give the functional form of the 
observable amplitude: 
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[16] Ay « Po « mGfo, 


where the ‘illumination factor’ G is the geometric mean of Gr and Gr, and 


fo = — Sin Bcos(y+e)|D'—hrcot'y| _ 
°" (D?-+h'cot?y)r*? (1—cos’B D?/r’)* ° 





[17] 


These relations apply, in the first instance, only to trails passing through S. 
Corresponding relations could be derived for other trails through the zone of 
observation, and they would not differ much from these if the zone were 
sufficiently small. For present purposes, then, [16] and [17] will be adopted as 
representative for all potentially observable meteor trails distinguished by the 
pair 8, y. 

If observations are confined to trails giving observed amplitudes equal to 
or exceeding some critical value A,, then they will be confined to meteors with 
masses equal to or exceeding some critical mass 


18] m, « A,/G fo. 


It is known empirically (5) that the number of meteors in a statistically 
large sample having masses equal to or exceeding any particular mass m, is 
inversely proportional to m,. (The relationship does not hold for all meteor 
showers, but it will be adopted here for want of anything better.) Accordingly, 
it is possible to detect only a fraction of the potentially observable trails 
passing through any part of the zone, and that fraction is proportional to 
G fo/A-. 
CONTOURS FOR CONSTANT GAIN 


If variations in the illumination factor can be ignored, then the observed 
fraction of the potentially observable trails does not vary with position in the 
zone of observation. Instead, it depends only on the defining pair 8, y, and has 
the functional form of fo. Since the number of potentially observable trails 
has the functional form of “in [8], there is a net ‘observability’ factor 


[19] fo a Fy = sin’B cot ¥[(F*—h*/sin*y)'—D]|D*—h'cot*y| 
0 ° PP +h/ (B+ #)siv’ pl —cos’8 D?/r}'[D? +h’cot* y] 


with h given by [4] or [5] and 7 by [11]. In this, the approximation 





[20] cot(y+o) ~ cot ¥/[1+h/(# +H#)siny] 
has been adopted; no significant error results when 2D is in the relevant 
range. 


A contour chart of Fo is given in Fig. 5, for a transmission axis of length 
2D = 1000 km., assuming “ = 100 km. as a representative meteor height 
and taking the earth’s radius # as 6370 km. This chart gives relative values, 
normalized to a possible maximum of 100, for the number of shower meteors 
counted per unit time when the radiant is in the 8, ¥ position, if the true flux 
density of meteors is constant. It will be seen that, as 8 and y change in the 
course of a day, there is ample scope for spurious observational changes in 
the numeric counts, beyond any random fluctuations in the true flux density. 
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Fic. 5. Normalized contours of the observability factor Fo, for isotropic antenna gain. 


A brief description of the zero contours may be useful. When |y| > 90°, 
Fy) = 0 because the shower radiant lies below the horizon and the trails are 
therefore formed in such positions that they cannot contain the center points 
of the corresponding Fresnel systems. At |y| = 90°, x = 90°; the zero in Fo 
results from the vanishing ionization density implied in the theory of trail 
formation, though this itself involves an inherent approximation which is 
invalid when x — 90°. When |y| < 2°, the radiant lies close to the vertical 
plane through 7°R, and the points of tangency lie well to the side of this plane; 
they lie, in fact, beyond the region of overlapping horizons, and hence the 
zero in Fo (represented by the hatched area in the figure). At |p| = 10.8°, 
Fy = 0 because cos y = 0 at the corresponding point S; this is not a true zero 
of Fo, since points P off the mid-plane will not lead to an exactly vanishing 
polarization factor, but a very low relative value of Fy is certainly to be 
expected. (The region of overlapping horizons is not very broad, at this dis- 
tance from the 7R axis, so cos y will approximate to zero throughout the 
relevant zone.) Finally, the decrease towards zero as |90°—8| — 90° results 
from two sin @ factors; one is introduced by the projection which gives the 
number of potentially observable trails, the other by the cos x factor in the 
ionization density. 


CONTOURS FOR A TYPICAL ANTENNA PATTERN 


When anisotropies in the antenna patterns are appreciable, variations in the 
illumination factor G must be taken into account. This could be partially 
accomplished by evaluating G for each S = S(y), and including it in [19] asa 
further factor. However, G may change appreciably within the relevant 
zone of observation, and it may actually limit the effective length of that zone. 
It therefore seems more appropriate to make an estimate of G as a function of 
distance from the mid-plane, at the height ”, and to integrate this value 
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along the length of the zone (at the distance h cot y to the side of the TR 
axis). This procedure would then provide the revised observability factor 


ad 
* 

21) Fo= J G dt -55--—S Fs op pss 
[21] : ao : r-+h/ (A+ )sinp][1—cos’8 D*/r*]'[D? +h'cot’y] 
with the length d given by [6] as before, and dz taken as an element of this 
length. 

Fig. 6 illustrates an illumination function such as might be established by 
five-element Yagis as transmitting and receiving antennae. Only the main 
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Fic. 6. Assumed antenna pattern construction lines, and contours of the corresponding 
illumination factor G, seen in projection on the horizontal plane through 7R. 


lobe of each is considered. It is taken to be beamed at the ” level above the 
mid-point, to fall to zero power on the horizon and at an equal angle above 
the beam center, and to fall to zero power at 50° to either side of the axis. 
Half power is attained at half the angles to zero power, and a normalized 
value of unity at the center is employed. The construction lines are shown 
on the left of the diagram, and refer to the same 2D, ”, and # values as 
before. The corresponding illumination factor is shown by contours on the 
right, and below that the projected positions of the observation zones (at 
y = hcoty) are indicated for several y’s. The increased imnortance of ohserva- 
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tion points P near the mid-plane is clear from this diagram. For comparison, 
the values of d, Gsd, and [G dz are tabulated after normalization to a maximum 
of unity (Table I). The contour chart of Fo* which corresponds to the assumed 
pattern is given in Fig. 7, after normalization to a maximum of 100. 


cos @ 








O68 





0-6 








0.2 











lvl 





Fic. 7. Normalized contours of the observability factor Fo*, for the illumination pattern of 
Fig. 6. 


TABLE I 


NORMALIZED EFFECTIVE LENGTHS OF THE ZONES OF OBSERVATION: d FOR ISOTROPIC ANTENNA 
GAIN; Gsd AND /Gdz FOR THE GAIN PATTERN OF FIG. 6 


7 10° 25° <20° .25°>30" 35° 40° 45° -S0° “60° -70° (80° 96° 


d 0.64 0.78 0.88 0.93 0.96 0.97 0.98 0.98 0.99 0.99 0.99 1.00 1.00 


1.00 
Gsd 0.00 0.04 0.15 0.26 0.36 0.47 0.57 0.65 0.72 0.78 0.88 0.94 0.98 1.00 
1.00 


JSGdz 0.00 0.01 0.18 0.36 0.53 0.64 0.71 0.77 0.82 0.87 0.92 0.96 0.99 


APPLICATION OF THE CONTOUR CHARTS 


Shower radiants move with the stars in the course of a day, and introduce 
rather involved time variations in the related values of the 6, Y coordinates. 
When observations are to be made on a number of showers over a particular 
path, the related theoretical predictions can be obtained most conveniently 
after a change to celestial polar coordinates. The radiant paths would then 
appear as circles on the corresponding contour charts, and these paths would 
be traced at a uniform angular rate. 

The transformation involves two parameters related to the mid-point of 
the TR axis, or to the point on the earth’s surface directly above it. These 
parameters are the latitude / of the point, and the azimuth a of the TR axis 
(measured to the east from north, say). For clarity, choose the azimuth either 
of T or of R, such that a lies between 0° and 180°; measure 8 from the direction 
so chosen, and then measure y > 0 in the hemisphere away from the pole. 
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The related quantities 


tan ¢) = tan a/sin I, 
[22] tan L = tan//sina, 
cos b = cosacos/ 


are useful subsidiaries, and are constants of the path employed. The relations 
cos P = cos8cosb + sin @sin db sin(L—y), 
[23] sin = sin 8 cos(L—y)/sin p 


may then be used to determine the polar distance p and the hour angle 
t = 7+% (after lower transit), corresponding to the coordinates 6, y. (See 
Fig. 8.) 
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Fic. 8. The celestial hemisphere over M, and the pertinent angles of the coordinate trans- 
formation. 


For some purposes—e.g., evaluating the effects of a finite shower radiant, 
or integrating to determine diurnal variations in random meteor counts— 
it is useful to plot the contour charts directly on a sphere. 
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ON THE DETERMINATION OF ELECTROMECHANICAL COUPLING 
COEFFICIENTS IN A POLYRESONANT PIEZOELECTRIC 
VIBRATOR! 


By E. A. G. SHAW 


ABSTRACT 
The well-known relationship between the electromechanical coupling coeff- 
cient and the separation (c) of resonance and antiresonance frequencies requires 
modification where several significant modes occur within a small frequency 
interval. An expression is given for the ‘‘true’’ values of o (i.e. those which 
would arise if each mode existed alone) in terms of the many moded system. 
The analysis neglects mechanical and electrical losses. 

In the course of a recent study of the resonant vibrations of prepolarized 
barium titanate ceramic disks, it became necessary to determine the coeffi- 
cients of electromechanical coupling in cases where several modes were 
separated by frequency intervals comparable with the differences between the 
frequencies of resonance and antiresonance of the individual modes. In such 
cases the characteristic variations of electrical impedance which occur in the 
resonance region are no longer related to the coupling coefficient in a simple 
manner. Indeed, in the most general case where there are many modes con- 
tained within a small frequency interval, each mode having a large coupling 
coefficient and large mechanical and electrical losses, the problem becomes 
highly intractable. The following treatment, however, while limited to loss- 
free systems, is capable of giving useful approximate values in many practical 
cases and has the merit of simplicity. 

By way of introduction, we consider first a piezoelectric vibrator having 
only one normal mode with appreciable electromechanical coupling. In stand- 
ard reference works (e.g. Cady (1)), it is shown that the electrical properties 
of such a resonator, neglecting losses, may be represented by the linear net- 
work of Fig. 1(a). (C; and Z, arise from electrical coupling with the mechan- 


Co 





(a) (b) 


Fic. 1. Equivalent electrical networks of a single mode piezoelectric vibrator: (a) exact 
form; (5) idealized form with impedance linear in ». 

t1Manuscript received May 19, 1955. 

Contribution from the Division of Applied Physics, National Research Council, Ottawa, Canada. 
Issued as N.R.C. No. 3686. 
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ical resonator, while Cy is the dielectric capacitance.) The network has zero 
impedance at the resonance frequency 


(1] v1 = 1/(2rVLiC ) 
and zero admittance at the antiresonance frequency 
[2] a, = nV 1+(Ci/C). 
Hence, 

[3] Ci/Co = (ay? — 4") /vy? & 20)/r1, 
where Oy = ay—1 KM. 


For the typical case of a longitudinally vibrating crystal, Mason (2) shows 
that the following relationship exists between o;/y; and the coupling coefficient 


ki: 
(4] ky? = = pilin fF 
4 V) 4 Vv) 


We now derive expressions for the impedances Z, and Z, of the two branches 
of the network 1(a) which have such a form that they may, with slight 
approximations, be made linear in »v. Combining [1] and [3] to eliminate C, 
and noting that Z) = (2rivCy)~', we have 

Zo = — (2miL,/v) (a1?—»,”) 


[5] ox —4niLyoy. 


Again, 
Z, = 2rivL,+(27ivC,) 


-1 
2rinLy (1 aan) = ¢ ae) | 
’ 1 


= txi(e—n)Li| 1 : =i cee | 


ll 


[6] ~ 4ri(v—v,) Li. 


The expressions [5] and [6] may be further simplified by expressing the ap- 
proximate forms of Zo and Z, in units of 4rL,0, = (2xv,Co)—'. Then we have 
Zo = —i and Z,' = i(v—v,)/o; as shown in Fig. 1(b). The real network of 
1(a) and the idealized network of 1() are approximately equivalent so long 
as v&v,; the divergences which occur at extremes of frequency will be 
considered later. 

The network of Fig. 1 may now be generalized as shown in Fig. 2, which 
represents a resonator having m resonances lying within a restricted frequency 
interval. It is obvious that the frequencies of zero impedance are ¥1,..., Ya, 
so that the resonance frequency of any one mode is not changed by the 
presence of the others. However, the frequencies of zero admittance a, ... , a 
are no longer equal to (v:+01),..., (%-+on), as would be the case if each 
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eects 


Cy Co 


cI a Be] 


Fic. 2. Equivalent networks of a polyresonant vibrator: (a) exact form; (b) idealized 
form. 


mode existed alone. Writing down the total admittance Y and setting it 
equal to zero, we have 


[7] Y= os ee +2) 








V-V, V—ve V—Vy 

_ I] (v—v,) toi (v—v2) (v—v3)...(V—Yn) +... + (V— 71) (v— v2)... (V—Yn-1) on 

oe ki eh ew, sy = = - = = nh = —_ 2 
I] (v—v,) 


= 0. 


Thus, ai,...,a@, are the roots of the numerator. Since the modes are num- 
bered in an arbitrary fashion, it will be sufficient to obtain a result in one 
case only, viz. for o;. To this end, all frequencies are referred to 11, the trans- 
formed values being distinguished with primes. Thus vp’ = y—y,, vy! = v1 —7, 
= 0, a,’ = a1—n, etc. Making these substitutions in the numerator of [7] 
and expanding with a view to obtaining an expression in descending powers of 
v’, we get 


Oe ii ant ui cae adit dn... 
=" — (v' 04’) (v’ 99"). . «(0 — Yai )on 
WET WW 5265 628g 8 +(—1)"vy/v0" 2. . vp’ 
—[or(’)” +. (HL) oe og... "J 
—[o2(v')” +. (1) yong’... 9] 


—[o,(v' P+. + (— 1) ae’... Vn-10n]. 
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This expression may be compared with the formal solution of [8]: 
n 
(9] I] @’-a,’) = ("+ ...+(-1)"ar'as’.. ah = 0. 
r=] 


It will be seen that, in both forms of the equation, the term in (v’)” has a 
coefficient of unity; hence the constant terms in [8] and [9] must also be 
identical. Furthermore, since v;’ = 0 by definition, only the second line of [8] 
makes any contribution to the constant term. Hence, 





/ / , / / / 
OW2V3 ~.-- Vn = AL Ag ..- Ay 
or 
ay/ats’ .. . On’ (ai1—1) (a@2— 1) . . . (&,— 11) 
[10] (= Ss a 
VoV3 «++ Vn (vo—v1) (vs3—1) . - «(Yn — 1) 
In the same manner it may be shown that 
pt ek el, ee ee ree pcsvees (emma 
” (v1 — Vm) (¥2— Vm) .. «(Ym—1— Ym) (¥m41— Pm) (¥m+2—Vm) “+ - (Yn — Ym) 


where m takes successive values from 1 to n. 

It will be noted that the right-hand side of [11] contains only observable 
quantities: the frequencies at which the network of 2(6) has zero impedance 
or zero admittance. Having thus obtained the values of ¢,, the corresponding 
coupling coefficients, k,, may then be deduced from [4]. However, the values 
so deduced are only accurate in so far as the network of 1(a) is accurately 
represented by 1(). A comparison of the impedance of Z;' of the idealized 
network 1(b) as a function of » with the corresponding value Z;/4rL,0, for 
the real network of 1(a) is made in Table I. It would appear that the diver- 











TABLE I 
Impedance 

v=0 vy = v,/2 y= 2y, y—> © 

Idealized network (Fig. 1()) —1iv;/o1 —tv;/20, +iv;/o1 +iv/o; 
Zi’ = (1/01) (v—”) 
Real network (Fig. 1(a)) —1@ —13v1/401 +13y;/4o; +iv/20; 
Bs ao (t/o1)( eh) 
4nL\o, — 2v 


gence in impedance between the two networks leads to computed values of 
om somewhat in excess of the exact values. Thus [11] provides excessive 
positive correction in respect of modes lying at frequencies very much greater 
than v,, and inadequate negative correction in respect of those at frequencies 
very much less than v». In most practical cases, however, only those modes 
lying close to vm are of much consequence, so that [11] gives o,, with sufficient 
accuracy for most purposes. 

It may be worth mentioning that values of ¢, calculated from [11] have 
been found to differ from a,—vm, by as much as a factor of 3 in the case of 
the group of modes lying in the thickness resonance region of barium titanate 


disks. 
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It would be exceedingly difficult to extend [11] to take account of mechan- 
ical and electrical energy losses, but it is possible to estimate their impor- 
tance by considering their effect on the electrical characteristics of a single 
mode resonator. Since this subject is discussed comprehensively in standard 
treatises (e.g. Cady (1)), only a few elementary results need be quoted here. 
While dielectric losses (those associated with Cy in Fig. 1(a)) may appre- 
ciably affect the magnitude of the antiresonance impedance, the frequency 
-of antiresonance is changed only to a negligible degree, even with such a lossy 
material as barium titanate ceramic. The mechanical losses, however, may 
modify appreciably both resonance and antiresonance frequencies, reducing 
the former and increasing the latter. Let a resistance Ry = 2mm,Lik; (ki = 
1/‘‘Q”’) be introduced in series with LZ; and C, in Fig. 1(a). Then it may be 
shown that, to the first order, 


[12] Z,/Zo| = kii/2o1, 


where Z, is the magnitude of the impedance minimum and Z, is the dielectric 
impedance alone. Furthermore, it may be shown that the difference o; between 
the frequencies of the impedance minimum and maximum is increased to the 
value o;* where 


[13] oi" = oi{1+ (Ry2n1?/2017) + eee ] =~ oi[{1+2|Z,/Zo|?+ see ]. 


Thus, provided the impedance minimum is a small fraction of the dielectric 
impedance, losses do not appreciably affect determination of the electro- 
mechanical coupling coefficient of a singly resonant vibrator. Likewise, it 
may be inferred that [11] will give an accurate value of c provided all of the 
modes lying close to ¢, have reasonably small impedance minima. 
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THE EFFECT OF A RANDOM NOISE BACKGROUND UPON THE 
DETECTION OF A SINUSOIDAL SIGNAL! 


By H. S. HEaps? 


ABSTRACT 


A sinusoidal signal is examined after reception upon a noise background of 
random phase and power distributed statistically according to a Rayleigh law. 
An expression is obtained for the probability distribution of the ratio of the 
power of the signal plus noise at any instant to the power of the noise which 
would be received at that instant in the absence of the signal. The corresponding 
result is given for the ratio of the averages over several observations. The equa- 
tions contain as a parameter the ratio of the signal power to the mean noise power. 
When each observed value of the power is the average over a small time interval 
the formulae are applicable provided the noise and the signal have the same 
frequency. A similar analysis is presented to deal with the case in. which the 
noise and the signal have different frequencies and in which each observation is 
the average over a small time interval. Comparison with the results of a previous 
paper, in which the signal was assumed to have a Rayleigh distribution in 
phase and power, indicates the effect of extreme fluctuation of an originally 
sinusoidal signal upon its resultant with a random noise background. 


1. INTRODUCTION 


A problem of frequent occurrence in connection with the propagation of 
acoustic or radio waves is the detection of an originally sinusoidal signal of 
the form asinwt after it has travelled through an inhomogeneous medium and 
has been added to unwanted interference during the process of detection. 
It is sometimes of interest to detect only a portion of the resultant signal 
such as, for example, the part reflected from a particular inhomogeneity or, 
in the case of radio waves, the part reflected from the ionosphere. 

If the receiver is situated sufficiently close to the source the signal upon 
reaching the receiver is essentially sinusoidal in character and has an envelope 
of constant amplitude a(¢) = @ with respect to time. As the receiver is placed 
at greater distances from the source the signal travels through more inhomo- 
geneities and its envelope a(¢) may differ significantly from that of a purely 
sinusoidal wave. For a receiver at a sufficient distance from the source the 
distortion of the signal may be such that its envelope contains no recognizable 
constant component and is best described statistically. Such a description may 
be necessary even if the inhomogeneities consist of quite small variations in 
the velocity of propagation and if the relative motions of the source, receiver, 
and elements of the medium are at low velocities. At any instant the sound 
reaching the receiver is the sum of waves which passed through the medium 
by paths of unequal travel time, and, if the travel times change by the order of 
a quarter wave period or more, a major change in the resultant wave will ensue 
(compare with Fig. 7 of Reference (12)). 

At the receiver the required signal is added to a noise background of dis- 
turbances generated within the receiver or the medium. This noise background 
may include unwanted components of the generated signal or of previous 

1Manuscript received April 21, 1956. 
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signals. Suppose that the noise background is slowly fluctuating in the sense 
that between two instants separated by a very small time interval the noise 
power remains unchanged while the phase increases by a periodic time factor. 
Suppose, however, that at two instants separated by a sufficiently large time 
interval the noise values may be related only statistically. At any instant both 
the signal and noise envelopes are specified by their amplitudes but the envelope 
of the signal plus the noise is dependent upon the amplitudes and the phase 
difference of the signal and the neise. 

Throughout the present paper the signal and the noise at any instant are 
represented in the complex form ptexp(i#) where p is real and is termed the 
power at that instant. It is assumed that at any instant all values of the phase 
6 are equally likely to occur. In the acoustical case it is known that noise 
due to reverberation from a large number of small random scatterers has a 
probability distribution of power according to a Rayleigh law (8, 13). Radio 
waves propagated through an inhomogeneous atmosphere and reflected from 
the ionosphere are found to contain a fluctuating component of power varia- 
tion closely approximated by a Rayleigh law (3, 7, 9). Also, in the detection 
of radio waves, a Gaussian distribution holds for the amplitude of fluctuation 
noise due to thermal agitation or to shot effect within the receiver (6). In 
view of these results, and by virtue of the central limit theorem that the 
superposition of a large number of random amplitudes has a probability 
distribution according to a Gaussian law regardless of the form of the individual 
amplitude distributions (11, p. 48), it appears reasonable to suppose that the 
assumption of a noise background of power distributed according to a Rayleigh 
law is valid in many instances. 

When the receiver is sufficiently distant from the source for the constant 
component of the received signal to be negligible the above considerations 
suggest that the power variation of the signal may also be supposed according 
to a Rayleigh law. The quantity detected by the receiver is then the sum of a 
Rayleigh distributed signal and a Rayleigh distributed noise. For this case the 
probability that, at any instant, the power of the received signal plus noise 
is a given multiple of the power of the noise alone has been obtained in an 
earlier paper (4), as has the corresponding result when the powers are averaged 
over several observations. The present paper considers the probability dis- 
tribution of the ratio of the power at any instant of a sinusoidal signal plus a 
Rayleigh distributed noise background to the power of the noise which would 
be received at that instant in the absence of the signal. The formulae obtained 
are thus applicable in the case of a signal whose envelope exhibits negligible 
fluctuation. Comparison of the results with those of the earlier paper indicates 
the effect of extreme fluctuation of the signal upon its impression on a random 
noise background. 


2. RATIO OF SIGNAL PLUS NOISE TO NOISE POWER FROM ONE OBSERVATION 


Kor a noise of phase and power distributed according to a Rayleigh law the 
probability, at any instant, that the noise power p, lies between p, and 
Patdp, with phase between @ and 6+dé6 in excess of the sinusoidal signal is 


[1] (1/2rp,) exp (— pn/ Pn) dp» 10 
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in terms of the mean noise power p,. The quantity dp,d@ represents twice the 
area p,'dp,'d@ swept by the end point of the amplitude vector p,%exp(i0) 
as p, and @ vary over the above ranges. 

When a signal of constant power p is added to a noise of instantaneous 
power ?, the condition for the noise plus the signal to be of power g is that @ 
satisfies the equation 


(2) cos 8 = (q—p—prn)/2pp,. 
This equation has two solutions in 6 if 
[3] la—p}| < pat < |gt+p? 


and no solution otherwise. As g varies between g and q+dq the end point of 
the vector p,'exp(i#), with @ determined by [2], ranges over an area of 


Thus the probability that the signal plus noise power lies between g and 
q+dq while the noise power lies between p, and p,+dp, is 
[5] (1/prn) exp (— pn/Pn) [4pq— (Pu—P—9)?] *dpndg. 


The probability that the signal plus noise power lies between g and g+dq 
is obtained by integration of [5] over the range of p, restricted by the condition 
[3]. The substitution p, = p+¢+2p}q! cos ¢ enables the integration to proceed 
over the range 0 < @ < z and leads to the expression 


(6] (1/pn) exp[— (P+9)/PnlIo(2p4q4/pr)dq 


in agreement with equations [3.10—-11] of Rice (11) and equation [5] of Blake 
(1) which refer to the corresponding distribution of amplitude. The expression 
[6] refers to the distribution of power and is plotted as a function of gq in 
Fig. 1 for p/p, values of 0.1, 1, and 2. The broken curves are of the corre- 


soll 


Pr 


2 


SUlln 





Fic. 1. Probability distribution P(g) of power g of the signal of mean power cp, plus 
noise of mean power pn. The continuous curves are for a sinusoidal signal (expression 6); the 
broken curves are for a Rayleigh signal (Reference 4). For each pair of curves the value of c 
is indicated. When c = 0.1 the two curves are indistinguishable. 
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sponding Rayleigh distributions for the case in which the signal is Rayleigh 
distributed of mean power p. The probability distribution of the phase of 
the resultant of a sinusoidal signal and a Rayleigh distributed noise has been 
investigated by Norton, Shultz, and Yarbrough (10). 

Putting g = kp, and dq = p,dk in [5] the probability that the power of 
the signal plus noise lies between k and k+dk times the power of the noise 
alone is P(k)dk where 


7] P(k) = (1/ xp) f exp(—Pe/a) MPP. (bubba) "TPP 
the integration being from p/(1+k?)? to p/(1—k?)?. The substitution p, = 


p(1+k—2kicos $)/(1—k)? reduces [7] to the form 
[8] P(e) = [p|1—k|-*/wp,] exp|— (1+) (1—k)*p/Pal 


x{ ate) i exp(a cos ¢)d¢@ — anf exp(a cos ¢) cos oae| 
0 0 


where a denotes 2k(1—k)-?/p, and |1—k| signifies the numerical value of 
1—k. The integrals appearing in [8] may be expressed in terms of the modified 
Bessel functions as rJo(a) and wJ;(a) respectively (14, p. 181). Substitution 
into [8], with c denoting p/j,, then leads to 
[9] P(b) = cll —kl-texp[— (142) (1—2)~*e]l (148) Z0(a) — 284 (2)]. 

For values of c « 1 the formula for P(k) takes the form 
[10] P(k) = c(1+k)|1—|-8 


as for a Rayleigh distributed signal when c < 1 (4). For values of c > 1 the 
asymptotic form of [9] bears no resemblance to the corresponding expression 





re 4 k6 ec 4c k & 


Fic. 2. Probability distribution of P(k), the probability that, at any instant, the signal 

lus noise power is equal to k times the noise power (equation 9). The dotted curves refer toa 
] ——— signal. On each curve the value of c, the mean signal power/mean noise power, is 
indicated. 

Fic. 3. Probability distribution of P(k) when c > 1. 
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for a Rayleigh distributed signal. When k = 1, a becomes infinite and the 
asymptotic formulae for Jo(a) and J;(a) (14, p. 203) reduce [9] to 
[11] P(1) = (1/8) (1/x)4c[exp(— 4c) ][1+-2/c]. 


The graphs of P(k) for c = 0.1, 1, 2, and >1 are drawn in Figs. 2 and 3, 
the dotted graphs being of the corresponding functions for the case of a 
Rayleigh distributed signal as calculated in Reference (4). The values of c 
are as indicated for each pair of curves. The most striking difference between 
the distributions for the sinusoidal and Rayleigh signals is a reduction of the 
maximum and, for c = 0.1 and 1, the occurrence of a minimum in ‘the sinu- 
soidal case. 

3. RATIO FROM THE MEAN OF SEVERAL RANDOM OBSERVATIONS 


Let the signal plus noise power be averaged over r+ 1 observations separated 
by sufficiently large time intervals for them to be independent although the 
mean characteristics of the signal and noise remain unchanged. Then the 
power of the signal plus noise averaged over the r+1 observations has a 
probability distribution dependent upon the value of r. Similarly the power 
of the noise alone averaged over the r+1 observations has a probability 
distribution dependent upon r. In this section the ‘mean’ noise power jn, 
refers to the average over a very long time while the ‘average’ noise power 
refers to the average of r+1 independent discrete samples. 

By repeated application of equation [5] it follows that the probability that 
the sum of the noise powers lies between p, and p,+dp, while the sum of the 
signal plus noise powers lies between g and g+dq is 


[12] dpadq( xia)" *exp(—Pa/Ba) f --- f bpm— (er p—9)"T? 


X [4p (2-91) — (x2—-x1 -— P—Y2t-y1)?]-4.... [40 (9-92) — (PnP —G+-97)71-4 
x dx dy, eee dx ,dy, 


where the integrations are over 0 < x1 <x2<...andO0O<y<y2<... 
subject to the condition that the expression within each square bracket of [12] 
remain positive. Repeated use of the identity 


t a 
13} ff t4mpy—(x—mp—yy I" [4m (ty) — (s—2—-p—t+y) | Mandy 
0 0 
= [ixrm™"*(m+1)™/(m—})p][4(m+1)pt—(s—mp—p—t)'T’, 
proved by application of the Faltung theorem for the double Laplace trans- 
form to the functions [4mpy— (x—mp—y)?]"-°/? and [4m(t—y) —(s—x—p— 
t+y)?]-? (see Appendix), enables [12] to be written in the form 
[14] (1/x)[r!/(2r) !](r+1)-"p-"B- [4 (7 +1) bg — (Pa -—1 —P— 9)". 
Putting g = kp, and dq = p,dk, the probability that the average signal 
plus noise power lies between k and k+dk times the average noise power is 
P,+41(k)dk, where 


[15] Pra(k) = (1/m)[r!/(2r) +1) "0h 


x f+ DepPA (bab P— hha) PaexD (— Pa/Pa)dPa 
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the integration being from (1+3)-?(r+1)p to (1—k4)-2(r+1)p. The sub- 
stitution p, = (1—k)-?(r+1)p(1+k—2k'cos ¢) transforms [15] to the form 


r+1-—r—1 


[16] Pr4a(k) = (1/m)[r!/(2r)!]2""k’ (r +1)" (1—k) pp, 


x exp[— (1+4)(1=#)+1)P/Pal| (1+) J expe cos ¢) sin’ ¢d¢ 


—2R' } exp(a cos ¢) cos ¢ sin’ as| 
0 


where a denotes 2k!(1—k)~-?(r +1) p/fn- The first integral in [16] has the value 
x} (4a)-’' (r+4)I,(a) (14, p. 79). The second integral may be obtained as the 
derivative of the first with respect to a. Substitution into the formula for 
P,41(k), with c denoting p/p,, leads to 

[17] P,(k) =k-» (1 —k)-*rcexp[— (1+k) (1—k)~?re][(1+k)I,-1(a) —2k4T,(a)] 


in which a is equal to 2k#(1—k)~*rc. 

For values of r between 0 and 11 and of a between 0 and 6 the values of 
I,(a) = (—1)"J,(ia) may be obtained from Reference (5), pp. 232-233. 
For large values of a the asymptotic form of J,(a) (14, p. 203) may be applied 
in the form 
[18] I,(a) = (2ra)-*(exp a)[1— (4r?—1)/8a+ (2r —1) (1+) (1—k})-2/2a] 


and leads to 

[19] P,(k) = 4(1/m)tk??-9)/4(1+k4)—*riclexp[—re(1+k})-?] 

X [1 — (4r?— 1) /8a+ (4r?— 1) (4r7?—9) /128a2+ (27 — 1) (1+k) (1—h#)-2/2a 
— (2r+1) (2r—1) (2r —3) (1 +k) (1—k?)-?/16a"]. 


For calculation of P,(k) when k = 1 the form [19] must be used and reduces 


to 
[20] P,(1) = (1/8) (1/r)*ric# [exp (— 4c) ][1+2(27—1)/rc]. 


The quantity P,(k)dk is the probability that the power of the sum of a 
sinusoidal signal and a Rayleigh distributed noise background averaged over r 
random observations lies between k and k+dk times the averaged noise power 
which would be obtained in the absence of the signal. The graphs of P,(k) 
for c = 0.1, 1, and >1 are drawn in Figs. 4-6 for r values of 1, 2, and 3. The 
probability that the ratio of the averaged powers exceeds k is given by 


(21] Q,(k) = J ”P, (k)dk. 


Numerical integration of the integral appearing in equation [21] leads to 
values of Q,(k) as drawn in Figs. 7 and 8. 

The median value of k with respect to P,(k) is the value of k at which 
Q,(k) is equal to 3. For a Rayleigh signal the median value of k is equal to 
c+1 regardless of the number of observations averaged (4), but for a sinusoidal 
signal it is apparent from Figs. 7 and 8 that the median value depends upon 
the magnitude of r and decreases with increasing r. For a given value of c 
the value of k at which Q,(k) is independent of r is less than the median 


value. 
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Fic. 4. Probability distribution of P,(&), the probability that the signal plus noise power 
averaged over r random samples is equal to k times the averaged noise power (equations 17, 
19). On each curve the value of 7 is indicated. c, the mean signal power/mean noise power, 


is 0.1. 
Fic. 5. Probability distribution of P,(k) when c = 1. The value of r is indicated on each 


curve. 
Fic. 6. Probability distribution of P,(k) when c >> 1. The value of r is indicated on each 


curve. 


Equation [9] of Section 2 compares the signal plus noise power with the noise 
power at a single instant and may be used to compare the averages over a 
small time interval provided the signal and noise have the same frequency. 
If the signal and the noise are of different frequency let wt denote the increase 
in phase difference between the noise and the signal during the time ¢ and 
suppose that the noise power does not change significantly throughout this 
time. An attempt to extend the approach of Section 2 to compare the signal 
and noise power averaged over the time ¢ with the power of the noise alone 
averaged over the time ¢ replaces the expression [4] by a more complicated 
function, and a similar subsequent analysis does not appear to be feasible. 
However there are two instances in which the resulting equations are much 
simplified. The first is when wt = 0, in which case the analysis proceeds as in 
Section 2 with p, p,, and g denoting averages over the time ¢; the second is 
when wt > 1, in which case the signal plus noise power averaged over the time 
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Fic. 7. Probability distribution of Q,(), the probability that the signal plus noise power 
averaged over r random samples is at least k times the averaged noise power (equation 21). 
The value of r is indicated on each curve. c denotes the signal power/mean noise power. 

Fic. 8. Probability distribution of Q,(k) when c >> 1. The value of r is indicated on each 
curve. 


t is given by Q = p+), where p and p, refer to averages over the interval ¢ 
(4). 

If wt >> 1 then p, = Q—p and the condition that the average noise plus 
signal power lies between k and k+dk times the average noise power is that 
pn lies between p, = p/(k—1) and p,+dp, where dp, = —pdk/(k—1)?. 
The probability of this occurrence is P(k)dk where 


[22] P(k) = c(k—1)~*exp[—c/(k—1)] 


and k is necessarily greater than unity. The probability that the average 
noise plus signal power is greater than k times the average noise power is 


[23] O(k) = J Pea sn ee 


If the above interval averages are again averaged over a random set of r 
time intervals the condition that the ratio of the average noise plus signal 
power to the average noise power lies between k and k+dk is that the sum of 
the noise powers lies between P, and P,+dP, where P, = rp/(k—1) and 


dP, = —rpdk/(k—1)?. Now the probability that the sum of r interval aver- 
ages of the noise power lies between P, and P,+dP, is 
[24] [1/(r—1)!]pn "Px 'exp(—Pn/Pn)dPn, 


and hence the probability that the ratio & lies between k and k+dk is P,(k)dk 
where 


[25] P,(k) = [1/(r—1)!](re)"(k—1)-?'exp[—rc/(k—1)] 
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Fic. 9. Probability distribution of P,() (equation 25) applying to the case wt >> 1, drawn 
for comparison with Figs. 4 and 5. The value of r is indicated on each curve. 

Fic. 10. Probability distribution of P,(k) when ¢ >> 1, drawn for comparison with Fig. 6. 

Fic. 11. Probability distribution of Q,() (equation 26) applying to the case wt >> 1, drawn 
for comparison with Fig. 7. a 

Fic. 12. Probability distribution of Q,(k) when c >> 1, drawn for comparison with Fig. 8. 


and & is at least unity. The probability that the average noise plus signal 
power exceeds the average noise power by at least a factor & is 


[26] @,(k) J "P, (k)dk 
1—expl—re/(-1)1D (1/n!)fre/(k—1)1", 
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which may also be expressed in the form 
Q,(k) = exp[—re/(k—1)] 20 (1/n!)[rc/(k—1)]”. 


Graphs of P,(k) and Q,(k) are drawn in Figs. 9-12 for c = 0.1, 1, and 
>l andr = 1, 2, and 3. It is to be noted that P,(k) and Q,(k) depend upon the 
values of r in contrast to the case of a Rayleigh distributed signal in which 
P,(k) and Q,(k) are independent of the number 7 of intervals considered 
(4). 

4. CONCLUDING REMARKS 


Comparison of the results of this and of the earlier paper (4) enables an 
evaluation to be made of the effect of fluctuation of the signal envelope upon 
its impression on a random noise background. Consider, for example, the case 
of a sinusoidal pulsed signal emitted by an underwater source of sound, it 
being assumed that the initial rise of each pulse occupies less time than is 
required for a significant fluctuation in the noise power. Suppose that a receiver 
placed near the source receives a signal of essentially constant envelope and of 
power large compared to that of the noise background. Suppose also that as 
the receiver recedes from the source the received signal power decreases in 
relation to the noise. 

Before arrival of the signal the receiver registers the noise background which, 
upon arrival of the signal, increases suddenly by a factor k. If the receiver is 
close to the source, the probability distribution of k is as shown in Fig. 6, 
and the effect of averaging the initial increase over several pulses is to con- 
centrate the distribution to the neighborhood of k = c. If, on the other hand, 
the signal plus noise is averaged over a complete pulse of duration ¢ such that 
wt >> 1, then the ratio k has the distribution of Fig. 10, which is similar to 
that of Fig. 6 except that k values of less than unity do not occur. 

If the receiver is sufficiently distant from the source for the mean noise 
background to have the power of the received signal, then the initial increase 
upon arrival of a signal has the distribution of Fig. 5. This differs significantly 
from the distribution of pulse averages shown in Fig. 9 in that the latter 
distributions indicate greater probability of k values between 1.5 and 2.5. 
For a detection system responsive to k values in excess of a threshold value 
between 1 and 2, Figs. 7 and 11 show that observation of pulse averages 
increases the probability of detection by about 10 to 20%. Comparison of 
Fig. 7 with Fig. 5 of Reference (4) shows that, with a threshold between 
1 and 2, the probability of detection of the initial increase is approximately 
the same as for a Rayleigh signal. But comparison of Fig. 11 with Fig. 8 of 
Reference (4) shows that, with a similar threshold, the probability of detection 
of a sinusoidal signal by pulse averages may be 30% more than of one that 
is Rayleigh distributed. 

If the mean noise power is 10 times the signal power, comparison of Figs. 
7 and 11 shows that a detection system with a threshold in excess of approxi- 
mately 1.2 detects fewer signals by pulse averages than by initial increases 
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caused by the arrival of a pulse. However, the detection probabilities are small 
in each case and differ by only a few per cent from the detection probabilities 
of a Rayleigh distributed signal. 
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APPENDIX 


The identity [13] may be proved as follows. Let J denote the expression 
t 8 
27) ff tampy—(e—mp—y) J" [4m (ty) — (6-2 pty)" aad. 
0 0 


By virtue of the Faltung theorem for the double Laplace transform J may also 
be written in the form 


27+ to 
{28] I = (2mi)~” SJ F,,(u, v) Fi (u, v) exp(us+vt)dudv 
Y— tao 
where 
[29] F,(u,v) = Jf tamoy— Ge—mp—y) "exp (—xu—yo)dxdy, 
The range of integration in [29] includes all positive values of x and y for which 
Ampy—(x—mp—y)? > 0. 
Putting x = mp+y—2(mpy)'cos @ equation [29] becomes 
[30] Fa(u,0) = exp(—mpu) f° (4mpy)"exp(—uy—oy)dy 


x f exp[(4mpy)*u cos 6] sin”” "ade 
0 


mT (m—})(4mp)”* (mpu’)*" [exp(—mpu)] 
x fo? Ina (Amputy)Sexp(—uy—ay)dy 
0 


by Reference (14), p. 79. The integral appearing in [30] is of a type listed in 
Reference (2), p. 197, and hence 

[31] Fm(u,v) = 24D (m—4)2?"-2(mp)"—!(u+v)-"exp[—mpuv/(u+v)]. 
Accordingly 


1+ to 


[32] I= (2mi)-*e*"P(m—3)2"™*(mpy" f exp|(t—s—mp—p)v]dv 
Y—to 


¥+ too 
x f “™" exp[(m+1) pv" /u] exp(su)du. 
Sy 


u 
lien 


The second integral in [32] has the value (2, p. 245) 
(2mi)[(m+1) pv?]-™/25™/2 Tm [4 (m+ 1) psv?]}, 
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and so [32] may be written as 
[33] I = (20i) a’? (m—4)2""-* (mp) 's"” 

¥+ too 

» 0" Im[4(m-+1)psv"]' exp[—4(m+1)psv’ exp (Tv) dv 

Y— too 
with T denoting t—s—mp—p+2[(m+1)ps]!. The last integral is equal to 
(2, p. 277) 
[34] [1/1 (m+ 4) ]x*[16 (m+ 1) ps}-™/?{ [16 (m+1)ps]}#7— T?} "4. 


Substitution of [34] into [33], and replacement of T by its value in terms of 
s and ft, finally leads to 


[35] I = [4mm™""(m+1)-™/(m— 4) p][4(m+ 1) pt— (s—mp—p—t)?]’, 
which completes the proof of equation [13]. 
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INTERNATIONAL COMPARISONS OF THE CANADIAN 
PRIMARY RADIUM STANDARD! 


By W. S. MIcHEL AND G. N. WHYTE 


ABSTRACT 


The content in milligrams of radium element of the Canadian primary radium 
standard was measured relative to the certified weight of a primary standard of 
the United States by an ionization method. The problem of correcting for dif- 
ferences in the attenuation of the gamma radiation in the body of each source is 
important, since the standards differ greatly in shape; a method of making this 
correction is described. Other intercomparisons involving the Canadian standard 
are discussed, and a mean value for the content of the standard is derived from 
owe comparisons. The agreement between this and the certified value is satis- 

actory. 


I. INTRODUCTION 


Measurements of the contents of radium sources, like measurements of 
mass and length, are based indirectly on a comparison with a national or 
international radium standard. In Canada all measurements of radium content, 
as well as the calibration of instruments for measuring gamma radiation, are 
based on the radium content of the Canadian primary radium standard, 
maintained in the laboratories of the National Research Council (and referred 
to hereafter as the NRC standard). An accurate knowledge of the value of 
this content is therefore of some importance. 

The NRC standard differs in two important respects from most other 
national standards”: the glass tube containing the radium salt is much smaller 
than those of the other standards, with the result that the salt fills most of the 
tube rather than a small part of it; and the certified value of the content is 
based not on a direct weighing of the salt but simply on a comparison with the 
former international primary standard on the basis of the ionization produced 
by the gamma radiation. It has been established (5) that no account was taken 
in this comparison of any difference in the wall and salt absorption in the two 
sources. Comparisons of the NRC standard with other standards that permit 
a calculation of the radium content of the former are thus of considerable 
interest. 

A welcome opportunity of this kind was provided by an interchange during 
the past few years of national radium standards for intercomparison purposes 
agreed on between the National Physical Laboratory (U.K.), the National 
Bureau of Standards (U.S.A.), and the National Research Council. Some of 
the results of these measurements have already been communicated in a joint 
statement by the three national laboratories to the Stockholm meeting of the 
Joint Commission on Standards, Units and Constants of Radioactivity. A 
joint statement on all the results obtained will be published elsewhere. 


1Manuscript received May 19, 1956. 
Contribution of the Division of Applied Physics, National Reseach Council, Ottawa, Ontario. 


Issued as N.R.C. No. 3697. 
*The history of the national and international radium standards has recently been outlined by 


Davenport, Mann, McCraven, and Smith (2). 
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This paper describes a recent comparison of the NRC standard and one of 
the standards of the National Bureau of Standards carried out in this labora- 
tory, and summarizes the results of a number of other comparisons involving 
the NRC standard. 


II. DETAILS OF THE NRC AND NBS STANDARDS 


The Canadian primary radium standard is a glass tube containing very 
pure anhydrous radium chloride, stated by the supplier (Union Miniére du 
Haut Katanga) to be free from thorium derivatives. It was acquired by the 
National Research Council in 1930 and subsequently certified by the Com- 
mission Internationale des Etalons de Radium to contain, in 1933, 24.24 mgm. 
of radium element. Standards measured at Paris and Vienna by the Com- 
mission are assigned Roman numerals: the NRC standard is number XIII. 

The United States primary radium standard (NBS standard) measured 
here is one of two acquired by the United States in 1936 from a set of 20 
prepared in 1934 by Professor Otto Hénigschmid. These standards were made 
up of carefully weighed amounts of highly purified anhydrous radium chloride. 
One of them is the present international standard in Paris, another is the 
present reserve standard in Vienna, and others have been acquired by the 
United Kingdom and Germany for use as national standards. The NBS 
standard used in this comparison was assigned the number 5440 by Hdnig- 
schmid. It was also compared with the former Paris and Vienna standards 
and received the number XV. 

Particulars of the two standards are given in Table I, and they are also 
shown in Fig. l(c) and (d). It should be noted that, while the two sources 


TABLE I 


PARTICULARS OF THE U.S. AND THE CANADIAN PRIMARY STANDARDS 











U.S. standard (NBS) Canadian standard (NRC) 





Identification or reference number 5440, XV XIII 

Radium content in mgm. element 20.45 2 24.24 }! 
20.38 3 

Internal diameter, mm. 34 1.5 § 

Wall thickness, mm. 0.27 0.25 

Length, mm. 38.3 10.5 § 


Nature of salt Dark brown grains Dark brown grains 





Certified value from the comparison with the primary standards of 1912 in Paris and Vienna, 


1933. 
2Calculated from the weighed value determined by Honigschmid on June 2, 1934; the relation 
1 gm. RaClz corresponds to 0.7612 gm. Ra was used. 

3Certified value obtained by comparison with the primary standards of 1912 in Paris and Vienna, 
corrected to 1934. 

44 more exact determination of the internal diameter of a portion of the tube was made as dis- 
cussed below. 

5As stated by Union Miniére du Haut Katanga. 


contain about the same amount of salt, the glass tube of the NBS standard is 
larger in diameter and longer than that of the NRC standard. If the standards 
are held vertically, the length of the salt column is about one-twentieth of the 
total length for the NBS source, as against four-fifths of the length for the 


NRC one. 
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Fic. 1. Photograph of four national radium standards. (a) U.S. standard number XIV 
(NBS); (6) U.K. standard (NPL); (c) Canadian standard (NRC); (d) U.S. standard number 
XV (NBS). (Photograph by courtesy of the Journal of Research of the National Bureau of 


Standards.) 
III. SELF-ABSORPTION CORRECTIONS 


The masses of radium element in radium preparations in equilibrium may 
be compared by comparing the gamma radiation emitted by them, since the 
intensity of this radiation is proportional to the mass of radium element 
present. Correction must be made, however, for differences in self-absorption— 
attenuation of the gamma radiation in the radium-bearing materia! itself 
(salt absorption) and in the walls of the container (wall absorption)—if the 
sources differ physically. Before making the ionization comparison, it was 
necessary to consider this self-absorption correction. 

Many of the old secondary standards are physically similar to each other 
and to the Hénigschmid standards. Even so, they may differ somewhat in 
grain size, in the relative amounts of free and occluded radon, and in grain 
position. Probably the differences in self-absorption due to these factors are 
minor; in any case, their effect has not been systematically investigated. The 
practice has been tu assume, at least in comparisons involving physically 
similar sources of approximately th® same content, that differences in self- 


absorption may be neglected. 
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When the assumption of equality of self-absorption cannot be made, the 
correction for each of the sources being compared must be known. Several 
workers have described their methods of making this correction (see, e.g., 
Perry (10) and Evans and Evans (3), who give numerous references). Since 
exact calculations cannot be made because of the difficulties in taking secondary 
effects into account, these methods are all empirical in part. In comparisons 
of standards it appears to be general practice to place sources in a horizontal 
position during measurement, care being taken to distribute the salt grains 
evenly along the bottom of the tube. However, it is difficult to calculate the 
correction for standards like the large ones in Fig. 1 in a horizontal position, 
since the salt consists of separated grains of uncertain size and position. 
Weiss (11) apparently made such corrections but gives little information on 
his methods of calculation. 

For the comparison of the physically very different sources described here, 
it seemed preferable to conduct the measurements with both standards held 
vertically. In this position a mean density for the salt can be calculated from 
a measurement of the volume it occupies at the bottom end of the tube, and 
a salt correction can be computed on the basis of this density. 

The salt and wall absorption corrections for the NRC standard held ver- 
tically, calculated from a knowledge of the absorption coefficients (9) valid 
for the National Research Council standard ionization chamber and using 
the value of 2.2 gm./cm.® determined at this laboratory for the density of the 
salt, are 0.66% and 0.23% respectively. The total self-absorption correction 
is therefore 0.89%. 

In order to obtain the corresponding correction for the NBS standard, one 
end of the standard was selected (all the ionization measurements were later 
made with this end down) and the salt shaken to the other end. The standard 
was then laid horizontal, immersed in Xylol, and viewed under a travelling 
microscope. Measurements of the internal diameter and wall thickness of the 
glass tube were made at intervals of 0.1 mm., starting from one end and 
working up to approximately the height that would be occupied by the salt 
in the vertical position (2.6 mm.). The wall correction for each zone was 
calculated using these measurements and the NRC linear absorption coefficient 
for soft glass of 8X10-* mm~. A mean wall correction was then found by 
weighting the correction for each zone according to the volume of salt in the 
zone; the result was 0.21%. In a similar way the weighted mean salt absorption 
correction was found to be 1.23%, giving a total self-absorption correction of 
1.44%, with an estimated error not exceeding 0.15%. This error is due mainly 
to the effects of refraction on the measurement of the inside dimensions of 
the tube. Incidentally, the density of the salt was found to be 2.2 gm./cm.', 
in exact agreement with that found for the NRC standard. 


IV. CONTENT MEASUREMENT 


‘The gamma-ray intensities from the two standards were compared on the 
NRC standard radium-measuring apparatus (9), which consists of a parallel- 
plate one-atmosphere ionization chamber with } in. lead walls lined with 
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aluminum, and a modified Townsend balance circuit for measuring the 
ionization current. During measurement, each source was held vertically on a 
strip of thin cellulose adhesive tape mounted in a light aluminum frame. The 
frame fitted into a holder placed 60 to 80 cm. from the ionization chamber. 
The holder could be moved by means of a vernier adjustment so that the axes 
of the sources were at the same distance from the chamber during the determi- 
nation of any one ratio. 

A total of 10 measurements, at different distances, of the ionization ratio 
of the two sources were made by three observers. In half of the measurements 
one standard was measured first, in half the other. Due precautions were taken 
to compensate for possible drift of the apparatus. The procedure also included 
rotation of each standard through 180° between successive readings. 

The mean gamma-ray ratio resulting from the measurements is 


NRC/NBS = 1.192+0.0019, 


where the uncertainty is the standard deviation of the mean of the measured 
ratios. Applying the absorption corrections calculated above, the content 


ratio becomes 
NRC/NBS = 1.186+0.0024. 


The weighed content of the NBS standard, corrected to 1933, is 20.46 mgm. 
of radium element. Based on this value, the content of the NRC standard is 


20.46 X 1.186 = 24.27+0.08 mgm.. 


where the estimated error includes an uncertainty of 0.03 mgm. in the weighed 
value of the NBS standard, based on the statement by Hénigschmid (4) that 
the total error in the weighed value of the sa/t could amount to 0.04 mgm. 


V. OTHER COMPARISONS INVOLVING THE NRC STANDARD 


Since its preparation in 1930, the NRC standard has been compared with 
other standards on a number of occasions. The results of several such com- 
parisons will be summarized briefly. 

(1) Direct Comparison with the Former Paris Primary Radium Standard 

The original certification of the NRC standard in 1933 was made by a 
direct comparison with the former (1911) international primary standard in 
Paris, the sources being laid horizontally on top of a flat ionization chamber 
not unlike the National Research Council chamber. (A comparison was also 
made with the ‘‘1911” reserve standard in Vienna, but details of the measure- 
ment are not available.) As described by Mme. Curie (1), the Paris standard 
consisted of 16.75 mgm. of radium element in the form of the chloride filling 
a length of 32 mm. in a glass tube 0.9 mm. in inside diameter and 0.27 mm. in 
wall thickness. On the basis of these dimensions and the assumption that the 
absorption coefficients for the Paris chamber do not differ much from those for 
the NRC chamber, one can calculate values of 0.20% for the salt absorption 
and 0.24% for the wall absorption, a total of 0.44%. The corresponding 
corrections for the NRC standard, taking into account the fact that in the 
horizontal position the salt does not occupy the full cross-section of the tube, 
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are 0.56% and 0.23%, giving a total of 0.79%. The content of 24.24 mgm. 
originally deduced for the NRC primary neglecting self-absorption corrections 
should thus be increased by 0.79—0.44 = 0.35%. The result is 24.32 mgm. 


(2) Indirect Comparison with the United Kingdom Primary Radium Standard 

During 1951 a glass source belonging to NRC and containing about 25 mgm. 
of radium was measured at both the National Physical Laboratory (NPL) in 
the United Kingdom and the National Research Council (6). Since the content 
deduced at NPL was based on the content of the U.K. primary radium 
standard (NPL standard), and that deduced at NRC was based on the NRC 
standard, the measurements constituted an indirect comparison of these 
standards. The fact that the particular glass source used was similar in form 
to the NPL standard simplified the comparison with the latter, since the self- 
absorption in the two could be assumed to be the same. The value obtained 
at the National Physical Laboratory was 24.02 mgm. The value obtained at 
NRC, recalculated on the basis of more recent absorption coefficients, is 
23.99 mgm. Taking the weighed content of the NPL standard to be correct, 
one may deduce that the certificate value of 24.24 mgm. for the NRC standard 
should be multiplied by a factor of 24.02/23.99, giving a value of 24.27 mgm. 
(3) Indirect Comparison with NBS Primary Standard 

In 1953 a platinum-enclosed secondary standard belonging to the National 
Bureau of Standards was compared at NRC with both the NBS and NRC 
standards (8). As in the previous case, this constituted an indirect comparison 
of the two standards. The absorption corrections were calculated in the 
manner described in Section III. The resulting content for the NRC standard 
was 24.17 mgm. 


(4) Direct Comparisons at the National Bureau of Standards 

Recently the primary standards of the United Kingdom, Canada, and the 
United States were intercompared by a number of methods at the National 
Bureau of Standards (2). Of the gamma-ray comparisons, that using the NBS 
standard electroscope lends itself most readily to analysis, since the same 
absorption corrections can be used as for the NRC chamber. 

The NRC standard was measured against the NPL standard and the two 
NBS standards, the sources being held in a horizontal position to one side of 
the chamber. By applying self-absorption corrections to the ionization ratios 
obtained, one can calculate the content of the NRC standard on the basis of 
the weighed contents of the others. Such corrections were calculated by 
assuming the salt to be fine-grained and to lie uniformly along the tubes, 
estimating the volume occupied on the basis of a density of 2.2 gm./cm.’, and 
determining the average absorption by graphical integration over the resulting 
segment-shaped cross-section. The assumption that the salt is fine-grained 
and uniform will certainly introduce some error, since the actual grains vary 
in size up to a few millimeters across and lie loosely along the bottom of the 
tube; but if one assumes the most extreme case—a single row of spherical 











MICHEL AND WHYTE: RADIUM STANDARD COMPARISONS 527 


grains along the length of the tube—the resulting salt correction is only lower 
by enough to make a difference of less than 0.2% in the content. The correction 
calculated for this extreme condition is certainly low, so that it seems unlikely 
that the corrections calculated on the assumption of fine-grained uniform salt 
are in error by much more than 0.1%. 

The content figures for the NRC standard deduced from the NBS measure- 
ments are presented in Table II, along with the other values obtained above. 
The mean value from all the ionization comparisons is 24.27 mgm. of radium 
element, with a standard deviation of about 0.1%. 


TABLE II 


VALUES OF THE RADIUM CONTENT OF THE NRC STANDARD 
DERIVED FROM IONIZATION COMPARISONS (AS OF 1933) 














Compared with At Year Derived content (mgm.) 

International primary Paris 1933 24.32 
NPL standard (5432, indirect) NPL & NRC 1951 24.27 
NBS standard (XV, indirect) NRC 1953 24.17 
NBS standard (XV) NRC 1953 24.27 
NBS (XIV) NBS 1954 24.29 
NPL (5432) NBS 1954 24.22 
NBS (XV) NBS 1954 24.33 

Mean = 24.27 





(5) Calorimetric Comparison at the National Bureau of Standards 

There exists one further measurement on the NRC standard. In the series 
of measurements at the National Bureau of Standards (2), the NRC source 
was compared with the NPL and NBS standards by a calorimetric method (7). 
The value derived from this work is 24.37 mgm., as of 1933. 


VI. DISCUSSION 


The mean value of 24.27 mgm. found for the radium content of the NRC 
standard from the ionization results in Table II differs by a little over 0.1% 
from the certified value of 24.24 mgm. Since the stated accuracy of the certified 
value is 0.2%, this agreement may be considered satisfactory. 

The calorimetric value of 24.37 is 0.4% higher than the value obtained from 
ionization measurements, and although it depends on fairly large corrections 
for the growth of polonium and radium E in the sources compared, the dis- 
crepancy is not understood. 

In any case, it is apparent that when one compares the NRC standard with 
a weighed standard by the ionization method and corrects for self-absorption 
with the aid of the NRC coefficients, one obtains a value for the content of 
the NRC standard that does not differ from the certificate value by more than 
about 0.3%. Conversely, if one deduces the content of another source from the 
certificate value of the NRC standard by this method, the result should be 
accurate to about 0.3%. This is the way in which the standard is used. 
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A METHOD FOR THE EVALUATION OF SOME LATTICE SUMS 
OCCURRING IN CALCULATIONS OF PHYSICAL 
PROPERTIES OF CRYSTALS. IP 


By G. C. BENSON AND H. P. SCHREIBER? 


ABSTRACT 


The methods of a previous paper are applied to the evaluation of lattice sums 
of the form 


bay _ krl*(k?+-22-+ m?)-* (r and s positive integers). 
kel m 


Useful transformations result when either 7 or s is even but not in general when 
both r and s are odd. A table of the function 


Si2) = te ki(k?-+2)-#/2 


for integral values of ¢ from 5 to 30i is sabiiaadl This permits the extension of the 
method to the casesr = 1,5 = landr=1,s =3. 


INTRODUCTION 


In a previous paper (4) (henceforth referred to as I) a procedure was de- 
scribed by which certain types of lattice sums can be transformed into fairly 
rapidly converging series of known tabulated functions (exponential, modified 
Bessel, and Riemann ¢-functions, etc.). The application of this method to 
triple sums of the form 


[1] T, (r,s) = pay RT (+2 +m’)™, 


where 7 and s are positive integers, is considered in the present paper. Sums of 
this type occur in calculations of elastic properties (2) (r = 2, s = 0; r = 2, 
= 2;r = 4, s = 0) and edge energies (3) (r = 1, s = 1) of crystals. 
It is shown in the next section that when 7 and s are both odd, application 
of the methods of I to T,,(r, s) gives rise to double sums which are not available 
as tabulated functions. The calculation of one of these, defined by 


[2] Silo) = Dike +P, 


is described in the third section and a table of values is presented there. 


TREATMENT OF THE SUMS 1T,(r, s) 


It will be assumed that the values of r, s, and m in equation [1] are such that 
T,(r, s) is a convergent sum (i.e. 2n > r+s+3). Application of the trans- 
formation described in I leads to the result 


a Ee — 4) Salts) 





[3] T,(r,s) = 


n—} 
+4 Kl (tm Ke-s2emJP FP) | 


kim 
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where, in general, 
[4] Si(r,s) = Lev (R +P)’. 
k,l 
Apart from S,_4(r, 5), the terms of equation [3] are all similar to those en- 
countered in the examples discussed in I, and the evaluation of 7,(7, s) by 


means of this equation is dependent upon the possibility of evaluating these 


new double sums. 
The functions S,(r, s) defined by equation [4] satisfy the symmetry relation 


[5] (7; 5) = SAS; r) ° 
The power r can be lowered by the reduction formula 
[6] S,(r,s) = S,-1(r—2, s) —S,(r—2, s +2), 
which after u successive applications gives the relation 
[7] Si(r,s) = Ud (—1)' "Cy Sines(r—2u, 5+2)). 
d 
In this equation, the quantities “C; are the usual combinatorial coefficients. 
For a further examination of S,(r, s) it is convenient to consider the fol- 
lowing possibilities: 
(i) one or both powers (7, s) even integers, 
(ii) both powers (r, s) odd integers. 


Case (1) 

When at least one of the powers is even, no generality is lost by putting 
r = 2p, where ? is an integer. This procedure is justified by the symmetry 
relation equation [5]. From equation [7] with » = p it can be seen that 


[8] S,(2p, s) = U8 (-1)? "Cy Sy-44(0, +24). 
i 
This equation expresses S,(2p, s) as a sum of terms of the type 


[9] SO, g) = Dalh@e+Py = Dire +r), 
where g is an integer. These, when rewritten as 
10 0.) = 4 ret eetee— pre |, 
can be treated by the procedures of I. 
It may also be noted that equation [8] with s = 0 gives 

(11] So(2p, OL (—1P] = Lo" (—1)? °C, Sy-v4 (0, 2) 
which leads to a number of useful relations such as 
[12] S,(0, 2) = $S,-1(0, 0), 

S,(0, 6) = 3[S,-3(0, 0) —3S,_2(0, 2)+3.S,_1(0, 4)] 
3[6S,1(0, 4) —S,-s(0, 0)] . 


ll 


[13] 








BENSON AND SCHREIBER: LATTICE SUMS. II 531 


Case (iz) 
When both the powers 7 and s are odd, the reduction formula gives 
[14] S(2p+1, s) = Lo (—1)? 7Cy Spss(1, s+2)). 
d - 


Terms of the form S,(1, 2g+1) satisfy the relation 
[15] S4(1, 29-1) (= 1) = Do (= 1)? Cy Spats, +1), 
J 


which is analogous to equation [11]. However they cannot be evaluated by 
application of the transformations described in I as can the sums S,(0, q). 

In connection with some calculations of edge energies of alkali halide 
crystals (6), it was necessary to compute values of 
[16] Sill, 1) = Fi(@) 
for a number of integral values of t. The method of computation and a table 
of results are given in the next section. Assuming that S,(1, 1) is known as a 
function of 7, S,(1, 3) can be obtained from 
[17] S,(1, 3) = 35,.(1, 1) 
but further extension of this process requires the evaluation of S,(1, 5), 
S,(1, 9), etc. Since there are no problems of physical significance giving rise to 
these higher sums, their evaluation does not appear to be worth while at this 
time. 

EVALUATION OF THE SUM S (@) 


The function defined by equation [2] converges for t > 4. However for low 
values of ¢ this convergence is very slow and it has not been possible to find a 
way of improving it. 

A simple, direct method for estimating lattice sums of this type has been 
described by Misra (5) and Born and Huang (2). Essentially it depends on 
splitting the sum into two parts. The first, composed of contributions from 
lattice points lying within a certain radius R, is evaluated by direct summation. 
Contributions from points outside R are approximated by integration. This 
procedure was adopted in the present case but a square array defined by 


k=1,2,..., Land/ = 1,2,..., ZL was found to be more suitable than a 
circular one. The sum over the points of this array is denoted by 
[18] SAL) = Dik +P)- 

kt 


An estimate of the contributions from the remaining points is given by the 


integral 


r=L+4 ro 
Sity= Jo ST aainnasiyt JS eye 
2 1 
[19] eke oe 


The required sum -“;,() is obtained from the equation 


[20] File) = lan [A (L)+4.(L)]. 
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The lowest value of ¢ considered was 5. Table I shows the results obtained 
for -4;(L) and 4,(L) at a number of values of L. Up to L = 20, the sums 


TABLE I 


SUMMARY OF CALCULATIONS OF S(e) 





i oR 


L No. of terms SL) JF (L) SAL) +4 s(L) 
in sum 
10 100 0.3954144222 0.0409723649 0.4363867871 
15 225 0.4085915330 0.0277817946 0.4363733276 
20 400 0.4153568757 0.0210129854 0.4363698611 
50 2500 0.4278340866 0.0085333016 0.4363673882 
60 3600 0.4292443118 0.007 1230023 0.4363673141 
70 4900 0.4302545430 0.0061127325 0.4363672755 
80 6400 0.4310138183 0.0053534354 0.4363672537 
90 8100 0.4316053150 0.0047619252 0.4363672402 
100 10000 0.4320791108 0.0042881209 0.4363672317 
© 0.436367 18* 


*Extrapolated value. 


S(L) were evaluated with a desk calculating machine; above L = 20 it was 
necessary to have the sum programmed and evaluated on an electronic com- 
puter. It is interesting to note that the sum -“;(L) is still changing in the third 
significant figure after adding 10‘ terms. The extrapolation of “.(L)+ 5(L) 
to L = was accomplished graphically by plotting log. {[-7s(L) +4 s(L)]/C} 
against the reciprocal of LZ. In this C is a constant approximately equal to 
S (0), 

Other values of .,() for integral ¢ from 6 to 30 were calculated in a similar 
way but the limit Z was reduced as ¢ increased, since the convergence of the 
sum became more rapid. The final results are summarized in Table II. The 


TABLE II 
VALUES OF THE sUM .,() = S- bl(k?-+/2)~ #2 
kyl 





t Silo) t S (2) 
5 4.3636718 X 10-1 18 1.9552101 x 10-3 
6 1.9249534 X 107! 19 1.3819966 X 10-3 
7 1.1162796 X 10-1 20 9.7697653 X 10-4 
8 7.1421138 X 10-2 21 6.9071868 X 10-4 
9 4.7814603 X 10-2 22 4.8836371 X 10-4 
10 3.2766933 X 10-2 23 3.4530381 X 10- 
11 2.2745505 X 10-2 24 2.4415708 X 10-4 
12 1.5905826 x 10-2 25 1.7264084 X 10-4 
13 1.1171254 X10? 26 1.2207360 x 10- 
14 7.8664535 X 107-3 27 8.6318215 X 10-5 
15 5.5480988 x 10-8 28 6.1035813 X 10-5 
16 3.9168060 x 10-3 29 4.3158667 X 10-5 
17 2.7668234 X 10-8 30 3.0517710 X 10-5 





graphical extrapolation was only required for the case ¢ = 5. It was also 
possible eventually to dispense with the integration step; thus results for 
t > 15 are based on the direct summation of 100 terms (i.e. L = 10). 
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In the case of ¢ = 5, the uncertainty in the graphical extrapolation is of 
the order of one unit in the eighth significant figure. For the sake of uniformity 
the results for ¢ > 6 were rounded to eight significant figures with a possible 
error of +5 in the ninth digit. It is of interest that for high values of ¢ the 
quantity log.-“;,(@) is a nearly linear function of ¢ with a limiting slope of 
—log.+/2. This is useful in interpolating at non-integral values of ¢. 


CONCLUSION 


In the previous sections it has been shown that 7,(r, s) can be reduced by 
the method of I to sums of known tabulated functions in all cases in which at 
least one of the powers r and s is even. When = 1,s = landr=1,s =3 
the sum can also be handled by this method in conjunction with the table of 
the function .7,(). In other cases where both r and s are odd, progress is 
hampered by the occurrence of certain double sums which cannot be reduced 
further and which have not been tabulated. It is easy to show that the above 
remarks also apply to sums similar to T,(r, s) but with k+/+™m restricted to 
odd or to even values. 

Other sums similar to T,(r, s) but with the limits for k and / extending from 
— © to +o instead of from 1 to ~ have been treated by Misra (5) and by 
Born and Bradburn (1). In this case the sum vanishes if either 7 or s is odd, thus 


[21] RNR +P +m)" =0  (rorsodd). 
k,i,m 
If both 7 and s are even the sum can be rewritten as 
[22] p Rl (k?+1?+-m’)™ = 4T,(r, s) (r and s even) 
k, i,m 


and the right-hand side evaluated as already described. 
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EDGE ENERGIES OF ALKALI HALIDE CRYSTALS! 


By H. P. SCHREIBER? AND G. C. BENSON 


ABSTRACT 


Theoretical calculations of the energies of <100> edges of sodium and po- 
tassium halide crystals are described. The present computations differ from 
those of previous workers by the inclusion of terms arising from van der Waals 
forces. It appears that these forces make a significant contribution to the edge 


energy. 
I. INTRODUCTION 


The heat of solution of a solid crystalline substance in a suitable solvent is 
a function of the particle size of the material. This dependence has been 
studied by several groups of workers (1, 6,8) and used as an experimental 
method for determining the surface energy of solids. Even in the ideal case of 
perfectly shaped crystalline particles of a uniform size, the difference between 
the results for fine and coarse material measured in these experiments should 
not be exactly proportional to the difference in specific areas of the samples 
since it must contain other contributions which, if the crystals are not too 
small, can be associated with edges and corners. At present this variation from 
linearity is beyond the range of the experimental techniques employed but 
some estimate of its magnitude can be obtained from theoretical models. 

Calculations of the edge energies of the alkali halides were made a number of 
years ago by Born and Stern (3) and by Lennard-Jones and Taylor (7). In 
these, the work necessary to form an edge was computed, assuming that the 
forces between ions could be represented as the sum of electrostatic and pure 
repulsive terms. More recently, Shuttleworth (13) has shown that van der 
Waals forces make a substantial contribution to the surface energies of the 
alkali halides. Although it is to be expected that this will also be true for the 
edge energy, verification of this point is of some interest. 

In the next section the main features of the theoretical model for the cal- 
culation of the edge energies are summarized. The necessary lattice sums are 
evaluated in Section III, and the results for the <100> edge energies of some 
alkali halides are presented in Section IV. 

II. DERIVATION OF THE EDGE ENERGY EXPRESSION 

The usual classical Born—Mayer model (12) of an ionic crystal is assumed. 
At 0° K. and neglecting zero point motion, the crystal is represented as an 
array of positive and negative ions fixed at the sites of a space lattice and 
interacting according to additive central force laws. 

In Fig. 1 an infinite crystal of the sodium chloride type, having nearest 
neighbor distance do, is represented schematically. An orthogonal Cartesian 
co-ordinate system is chosen to coincide with the principal crystalline axes. 


1Manuscript received May 17, 1955. 
Contribution from the Division of Pure Chemistry, National Research Council, Ottawa, Canada. 
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Fic. 1. The co-ordinate system. Diagram shows the XOY plane of an infinite sodium 
chloride type crystal with nearest neighbor distance ao. The dashed lines indicate division 
of the crystal by the planes x = —}a9 and y = —}apo to form four parts. 


When the crystal is divided by the planes x = —}a9 and y = —3a, (broken 
lines in Fig. 1) and the four parts infinitely separated, four <100> edges 
(perpendicular to the page in Fig. 1) are created. 

If, following Born and Stern (3), the four parts are designated by the 
numbers 1, 2, 3, 4, the mean edge free energy per unit length is given by 
[1] k = (Uist U2) (4Z)7 
where U;; and U2, represent the energies of interaction between ions in diag- 
onally opposite parts and L is the length of one of the edges. Stranski (14) has 
discussed the definition of the edge energy in unsymmetrical cases where U3 
and Us, are not equal. However the formulation of Born and Stern is quite 
adequate in the present case, since U;3 and U2, are equal, and thus the energy 
associated with unit length of any one of the four edges is 
[2] « = Uj;(2L)—. 

From a consideration of the periodicity occurring along the edge (i.e. in the 
Z-direction) it is evident that 
[3] Uy3L~ = (U;+U_)(2a9)! 


where U, and U_ represent the interaction energy of all the ions in part 1 
with those of part 3 which lie in the planes Z = 0 and Z = ap respectively. 
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The co-ordinates of an ion in part 1 are p10, pedo, P3do where fi, po, p3 are in- 
tegers such that 


ti > 0, Pp > 0, and —-w~ Kg pj<g om. 


An ion of part 3 lying in the plane Z = 0 has co-ordinates —qido, —q2do, 0 
where g, and q are integers satisfying the conditions gq; > 1, gz > 1. The 
distance between these two ions is 


[+] R(p1, pa Ps, Qs G2) = ol (Pitgi)?+(po+qe)?+ps?]}. 


The energy of interaction of two ions a distance R apart is represented by 
€,,(R), €_(R), or €,_(R) according to the respective charges of the inter- 
acting ions. Assuming the origin to be at a positive site (as shown in Fig. 1), 
5} Usa YF &R+ FY €M+ LF €,_(R) 


PitDatpseven Pitpatpsodd PitPatPsteitgsodd 
q@it+qseven @it¢s odd 


where the sums extend over all permitted values of the integers pi, po, p3, qi, 
and g2 consistent with the restriction listed below each summation sign. By 
shifting the origin to the negative ion at (0, 0, ao) it can be seen that 


6] U= DY G&R+ LY €H+ YF €,_(R) 
PitDatPs odd PitPatPs even PitDatPstaitds odd 
@i+@_ odd @itgs even 


and hence 


(7] c= aayf > (Gi (R) + E_(R)) 


2 
PitDatPsatgitds even 


g(r) | 


PitDatPstqites odd 
These sums can be simplified by the substitutions 
Attn = k, PotQ = i, p=m 


where the &, /, and m must be integers such that k>1,/> 1, and 
—«o < m< o. The interionic distance R is then written as 


[8] R(p, po, P3, qi, 2) = r(k, i m) = alk? +1?+ m?}}. 
Since any fixed value of k can be made up from & pairs of values of p; and q;’ 
i.e. (0, &), (1, R-1),..., (R—-1, 1) and similarly any fixed / from / pairs of 


values of p2 and q, it is apparent that the edge energy can be written as 


[9] «= 200)" 37E7 kl (G.4(r)+E__(r)) + D1 dow ae, | 
et bem am at lim odd 
Following Huggins and Mayer (5) the interaction €;;(7) between two ions 


of species 7 and j at distance r is taken to be 


2 
i d —f 
[10] Ciy(r) = 22,5 -Y-G ob bife 


This is the sum of Coulombic, van der Waals attractive, and exponential 
repulsive terms, in which Z; and Z; are the valencies of the ions and e the 
electronic charge. For alkali halide crystals the van der Waals coefficients c;; 
and d;; have been tabulated by Mayer (9), the repulsive constants 0, b;, and p 
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by Huggins and Mayer (5), and excluding lithium salts, the Pauling factors* 
S++, f—, f+- are 1.25, 0.75, and 1.00 respectively. 
After the aaa for the &;;(r) is substituted, equation [9] becomes 


{ll] «= ea) # Bis A(,_py ee) Ds") 


heen pp Siatbtes) Ds") 
On amp 2 Tey 2 
sof hii +O7m. exp( vise), 2b,b_ exp( zing) 


/ 





repulsive terms arising from the interactions of ions farther apart than third 
neighbors having been neglected. The symbols B, and D,’, etc., are defined by 
the equations 


[12] B, a *) yo (-1) el (1 4m’) 4, 
[13] ~ 2 yo RR +P +m’), 

abise odd 
[14] D2!" = dad 1 Dow RI(k?+1°+-m?)™. 


k+l+m even 


Although the sums D,,’ and D,,’’ are convergent only for m > 5/2, a uniform 
approach to the evaluation of these sums and of B, is obtained by expressing 
the latter in a formal way as 


{15] B, = D," —Dy'. 
III. EVALUATION OF THE SUMS 8B,, De’, De’, Ds’, AND Dy” 
Application of the transformations described by van der Hoff and Benson (4) 
and Benson and Schreiber (2) to the sums De,’ and D2,”’ leads to the results 


[16] Dd = Pes T'(n—3), Fon 1(@ ) 


n—} eee 


Bee odd 


anes m ™m ‘aie at | 
+ 2 da(-)) ce K,-4(4mr/k' +1) 
k+l even 


and 


[17] Dey” = 2 p(n 4) Fm-1() 





T'(n) 
@ ~ TT n-t 
+2 da wa om Ky-4(armrn/k° +l’) 
k+l even 


*These are usually denoted by c,,, c__, and cy_. The present notation has been adopted to avoid 
confusion with the van der Waals constants. 
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where the function 


[18] Silo) = Lyrae +P)” 
k,l 


has been tabulated for a number of integral values of ¢t in Reference (2). 
After subtracting equation [16] from equation [17] and putting m = 3, the 
equation 


-_ By = 40258 (=1) "el Ko(rm VE FP) 


k,l m odd 


is obtained. Essentially this same result has been given by Born and Stern (3) 
but their method of derivation differs from that outlined above. 

When x is an integer the modified Bessel functions in equations [16] and [17] 
are of half integral order and can be expressed in terms of exponentials. The 
formulae for D2,’ and D2,’’ then become 


> 
20] Day’ = = ee BN Suyi(@) 


'(n—1)!L2"(n—1)! ° 
ee (ss Litt 2a di (-"7,) | 
k+l odd Sehane 
and 
n ! 
[21] Ds = o- “Kaa -1)! | ge= i "on—1( © ) 
f Song riri(-»r.) | 
bee inen tea odd 
where 
kim n—1 (n—1+y)! 1 


22] To GP ex rmr/k +l") 20 Gi») Gam EEE) 


Numerical computation of the lattice sums B,, Dg’, De’, Ds’, and Ds” using 
these equations yields the values* 
B, = 0.021822., 
D,’ = 0.230288:, D,'’ = 0.285893, 
D;' = 0.0368379¢, Ds" = 0.0748789s. 


IV. RESULTS OF ee LATIONS FOR SODIUM AND POTASSIUM 
ALIDE CRYSTALS 


The contributions of the various terms of equation [11] to the <100> edge 
energies of the sodium and potassium halides are displayed in Table I. For 
consistency in these calculations, the electronic charge was taken to be 
4.77 X10-'° e.s.u. since this was the value employed in establishing the po- 
tential constants (see References 5 and 9); furthermore values of a) calculated 
from the potential constants by Shuttleworth (13) were used rather than 
experimental lattice spacings. 


*Less accurate values which have been published by other authors are: 
B, = 0.02186, Born and Stern (8); Ds’ = 0.0865 and D,'' = 0.0745, Lennard-Jones and 
Taylor (7). 








f 
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TABLE I 


CONTRIBUTION OF THE DIFFERENT FORCES TO THE EDGE ENERGY OF AN UNDISTORTED 
<100> EDGE OF SOME ALKALI HALIDE CRYSTALS 


7 Energies in 
10-° erg cm.7! NaF NaCl NaBr Nal KF KCl KBr KI 


Electrostatic 4.734 3.212 2.857 2.434 3.522 2.592 2.348 2.048 
Dipole—dipole —0.550 —0.756 —0.813 —0.895 —0.571 —0.594 —0.608 —0.644 
Dipole—quadrupole —0.027 —0.047 —0.052 —0.061 —0.019 —0.025 —0.027 —0.033 
Repulsive 0.750 = 0.441 0.409 0.388) = 0.343 0.140 0.123016 
Total x 4.91 2.85 2.40 1.87 3.27 2.11 1.84 1.49 
Born and Stern 5.08 3.42 3.05 2.60 3.82 2.75 2.50 2.18 
Lennard-Jones 
and Taylor 5.13 3.97 3.95 3.06 


V. CONCLUSION 

For comparison, edge energies calculated from the formula derived by 
Born and Stern* (3) and the results of Lennard-Jones and Taylor (7) are shown 
in the last two lines of Table I. Although different forms of repulsive potential 
have been used in the three calculations, the main difference’ between the 
present results and the other two sets can be attributed to the lowering effect 
of the van der Waals terms. The contributions from these forces to the total 
edge energy vary progressively from 12 to 50% for the sodium halides and 
from 18 to 45% for the potassium salts. In general, van der Waals forces play 
only a minor role in the lattice energy, but a significant contribution to the 
surface energy has been noted by Shuttleworth (13). The relative contributions 
of van der Waals forces to the lattice, surface, and edge energies of sodium 
chloride are illustrated in Table II. 


TABLE II 

COMPARISON OF THE RELATIVE CONTRIBUTIONS OF DIFFERENT 
FORCES TO THE LATTICE ENERGY, {100} SURFACE ENERGY, 

AND <.100> EDGE ENERGY OF SODIUM CHLORIDE 


{100} surface <100> edge 








Lattice energy, energy, energy, 
kcal. mol.~! erg cm.~? 10-° erg cm.7! 
Electrostatic 207.2 345.6 3.212 
Dipole-dipole 5.7 79.6 —0.756 
Dipole—quadrupole 0.7 8.6 —0.047 
Repulsive — 29.0 — 245.3 0.441 
Zero point — 1.7 
Total 183 189 2.85 
% Contribution of 
van der Waals terms 3.5 46 28 


*Born and Stern give two alternate formulae; one in terms of lattice parameter, the other in terms 
of density. The numerical factor in the latter is incorrect. The values of the edge energy which 
Lennard-Jones and Taylor attribute to Born and Stern are calculated from this erroneous equation. 
It should also be noted that in computing the last two lines of Table I, experimentally determined 
interionic distances were used, which differ slightly from Shuttleworth’s values of ao. The effect of 
this difference on the edge energy is slight. 
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Finally, it should be noted that along with previous authors it has been 
assumed inherently in the derivation of equation [11] that the ions of the 
newly formed edge are unpolarized and remain fixed at the regular lattice 
sites. The problem of calculating the equilibrium configuration of the ions 
near a crystal edge is a very difficult one and is complicated by the fact that 
as yet no really comprehensive treatment has been given for the distortion at 
the surface of the crystal. Some calculations by Moliére, Rathje, and 
Stranski (10) indicate a tendency for ion pair formation in the crystal surfaces 
and edges of some alkali halides. However Patterson (11) has shown that these 
configurations in the surface are unstable when the effect of van der Waals 
forces is considered. Since polarization and distortion effects undoubtedly 
occur in the vicinity of the edges of real crystals, the values in Table I must be 
considered as first approximations to the edge energy. 
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THE FISSION YIELDS OF ISOTOPES OF XENON AND KRYPTON 
IN THE NEUTRON FISSION OF U** AND U**! 


By R. K. WANLEssS AND H. G. THODE 


ABSTRACT 


Several uranium samples have been irradiated under various conditions in the 
Chalk River N.R.X. reactor and in the Los Alamos fast reactor. The fission 
gases, xenon and krypton, have been extracted from the irradiated material and 
the relative isotopic abundances determined mass spectrometrically. Fine 
structure in the mass — fission yield curve has been found in both the xenon and 
krypton mass ranges for U** + n and U2 + n fission. It is observed that this 
fine structure shifts to lower masses, as predicted, for neutron fission of U**, 
The proportions of U** and U** fission that have occurred in the irradiated 
samples have been estimated from the percentage change in fine structure. 


INTRODUCTION 


Early mass spectrometric determinations of the relative fission yields of 
the stable and long-lived isotopes of krypton and xenon (10, 13), produced 
in the thermal neutron fission of U**5, established the fact that the yields of 
Xe! and Xe" fell considerably above the smooth curve drawn through the 
adjacent xenon isotope yields. Glendenin (6) has postulated a chain branching 
mechanism to explain this ‘‘fine structure’ in the fission yield curve. He sug- 
gested that primary fission fragments formed with 51 and 83 neutrons (one 
neutron in excess of a closed shell) would have sufficient excitation to emit 
the loosely bound 51st or 83rd neutron rather than achieve stability by under- 
going 8 decay. This mechanism provided only a qualitative explanation of 
the observed yields in U™* + n fission. Wiles et al. (15, 16) suggested that 
better quantitative agreement between predicted and observed yields could 
be obtained if a primary fission preference for nuclei with 82 neutrons was 
assumed in addition to the Glendenin chain branching mechanism. 

Recently Pappas (11) has extended the ‘‘neutron boil-off’’ mechanism to 
include neutron emission by primary fission fragments with one, three, five, 
or seven neutrons in excess of a closed shell. By also taking into account a 
primary fission preference for 82-neutron nuclei, first suggested by Wiles, and 
tentatively assigning the 1.52-sec. delayed neutron activity to Sb, Pappas 
obtained a reasonably accurate quantitative explanation of the observed 
U*% + mn fission yields. 

The ‘‘neutron boil-off’’ mechanism, as modified by Pappas, can be used to 
predict fine structure effects for other fissioning nuclei. The fission yields of 
the xenon and krypton isotopes in neutron fission of U™* are reported in this 
paper and compared with the yields predicted by Pappas. 

In the earlier investigations of the relative yields of the xenon and krypton 
isotopes in U** + n fission, small variations in the abundances of the isotopes 
were often observed, depending on the irradiation conditions, position in 
nuclear reactor, etc., of the normal uranium samples. In particular, the high 


1Manuscript received May 9, 1958. c 
Contribution from Hamulton College, McMaster University, Hamilton, Ontario. 
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yield at mass 134 was found to vary by a significant amount, the lowest yield 
being obtained when the samples were irradiated in regions of higher neutron 
flux, and greater neutron energies. Since the energy distribution of pile neutrons 
was such that some U*88 + n fission occurred, it was felt that the variation in 
the 134 mass yield found from sample to sample was due primarily to the 
different proportions of U** and U** neutron fission which had occurred. 
Although the mass-yield curve for U** fission is known to parallel that for 
U*5 fission fairly well in this mass range, the fine structure or abnormal yield 
at mass 134 was not expected to be the same for the two nuclei. In order to 
investigate these variations further and to obtain a satisfactory comparison 
between the relative fission yields for U*® and U* fission, samples of natural 
uranium and samples of uranium depleted in U** were irradiated in the 
thermal column of the N.R.X. pile and in the Los Alamos fast reactor re- 
spectively. The xenon and krypton gases extracted from these samples gave 
the isotope fission patterns for the strictly thermal neutron fission of U*** and 
for the fission of U*88 (average neutron energy 0.5 Mev.). 

Also in the course of this study, samples of natural uranium were irradiated 
under widely different conditions. For these samples, it was possible to esti- 
mate the proportion of U* and U8 fission which had occurred by a com- 
parison of the xenon isotope fission pattern of the natural uranium with those 
obtained for U*** (thermal neutrons) and U?* fission alone. 

It is realized, of course, that the fission isotope patterns obtained for U** 
will to some extent depend on the energy distribution of the neutrons during 
irradiation since symmetric fission increases as the energy of the neutrons 
increases from thermal energies to 14 Mev. (12). Also, fission yields will depend 
on neutron flux where nuclides with high neutron capture cross sections are 
involved. Therefore, any variations in yields due to these factors will be 
superimposed on effects due to varying proportions of U*** and U* fission. 
Results indicate that although the varying proportions of U*** fission produce 
by far the greater effect on the fission yields in the xenon range, there is some 
indication that the distribution of fission products for U** varies significantly 
as we go from a thermal distribution of neutrons to neutrons having an average 
energy of 0.5 Mev. 


EXPERIMENTAL 


(i) Irradiation of Samples 

In order to more fully investigate the small variations in fission yields 
observed in the early work, old samples of fission product xenon and krypton 
were repurified. These samples, from uranium irradiated under widely differing 
conditions, were carefully compared by analyzing each sample on the mass 
spectrometer under identical operating conditions. 

The irradiation data for the new samples used in this work are summarized 
in Table I. Irradiation 1 was normal uranium metal irradiated in the thermal 
column of the N.R.X. reactor at Chalk River. In this position the neutron 
energy distribution should be entirely thermal and hence no U?* fission would 
be expected. Irradiation 2 was normal uranium metal sheathed in cadmium 
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TABLE I 
IRRADIATION DATA 





Approximate 





Irradiation Description Reactor Approximate integrated 
number neutron flux,* neutron flux, 
n./cm.?/sec. n./cm.? 
1 Normal uranium metal N.RUX,, 5.0 10" 2.910!7 
Chalk River (thermal column) 

2 Normal uranium metal N.R.X., 1.510" 1.51019 
sheathed in cadmium Chalk River 

3 Normal uranium metal Fast reactor, 1.010" 5.6 X 108 
Los Alamos 

4 UO; depleted in U* Fast reactor, 1.0X10" 1.710" 
Los Alamos 





* Flux values for the Chalk River irradiations are the estimated thermal neutron fluxes. Los Alamos 
irradiations were carried out in a fast neutron flux with average energy of 0.5 Mev. 
1/16 in. thick and irradiated in a high flux position in the N.R.X. reactor. 
Since cadmium absorbs most thermal neutrons and since there is an appreciable 
fraction of high energy neutrons at this irradiation position, considerable U**8 
fission would be expected in this sample. Irradiation 3 was normal uranium 
metal irradiated in the fast reactor at Los Alamos where the number of thermal 
neutrons is negligible. This sample should therefore show even more U** 
fission than sample 2. Irradiation 4 was UO; depleted 30 to 50 times in U** 
and was irradiated in the Los Alamos fast reactor. This sample should contain 
a negligible contribution of U** fission. Duplicate samples were irradiated in 
each case to provide checks on the extraction, purification, and mass spectrom- 
eter analysis of the rare gas fission products. 


(ii) Preparation of Rare Gas Samples 

The apparatus used in this laboratory for the extraction and purification 
of fission product gases from irradiated uranium metal has been described 
by Arrol, Chackett, and Epstein (1). The method consists of dissolving the 
uranium in a saturated solution of cupric potassium chloride, drying the 
released gases with solid potassium hydroxide and magnesium perchlorate, 
and finally, purifying with calcium vapor. This procedure was followed exactly 
for samples 1, 2, and 3. In the case of sample 4 a slight modification was re- 
quired since this sample was in the form of finely powdered UO;. The UO; 
was dissolved in concentrated sulphuric acid rather than the cupric potassium 
chloride solution used for the uranium metal samples. The subsequent handling 
of the released gases was the same as described in Reference (1). 


(iii) Mass Spectrometry 

All gas samples were analyzed on a 180° direction focusing mass spectrom- 
eter. Magnetic scanning was employed throughout to focus ions of 1500 v. 
energy. An Applied Physics Corporation vibrating reed electrometer was used 
to amplify the ion currents which were recorded on a Leeds and Northrup 
Speedomax strip chart recorder. Typical isotope patterns of fission product 
krypton and xenon are shown in Figs. 1 and 2 respectively. At least 10 double 
spectrograms were recorded and averaged for each analysis. 
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Fic. 1. Mass spectrogram of fission pro- Fic. 2. Mass spectrogram of fission pro- 
duct krypton. duct xenon. 


Care was taken to ensure that there were no residuals in the mass spectrom- 
eter before samples were admitted, and each sample was checked for im- 
purities in the neighborhood of the krypton and xenon mass ranges. To remove 
any condensable impurities that might have been present in the samples a 
liquid air trap was employed between the sample input line and the mass 
spectrometer. The isotope ratios obtained for the xenon and krypton isotopes 
were reproducible to better than +0.5%. 


RESULTS AND DISCUSSION 


(i) Xenon Isotope Fission Yields 

The mass spectrometer measurements of fission product isotope abundances 
only give relative fission yields for the isotopes of the element being studied. 
These yields must be normalized in some manner before they can be compared 
with yield data for other elements or with yields for the same element for a 
different fissioning nucleus. We have assigned the value 3.00%* to the fission 
yield of Xe! in U?* + m fission (thermal column irradiations). 

The fission yield data for the old fission product samples which were re- 
analyzed are given in Table II (normalized at 3.00% for Xe"). The approxi- 


*Recently reported independent determinations of the yield of the 181 mass chain are: 
Katcoff and Rubinson—2.97% based on an absolute determination of the yield of Xe (8) and 
the relative yields of Xe™ and Xe™ ane by Macnamara et al. (10), and Wiles et al. (16). 
Yaffe et al.—3.21% based on an absolute 
yields of Ba and I™ (2). 


termination of the yield of Ba‘ (17) and on the relative 
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TABLE II 
FISSION YIELDS OF XENON ISOTOPES UNDER DIFFERING IRRADIATION CONDITIONS 








Approximate % Fission yield 
Approximate integrated —- 
Sample neutron flux, neutron flux, 

n./cm.?/sec. n./cm.? Xe3! Xe! Xel4 Xels6* 
A 5.0 X 108 2.0 X 108 3.00t 4.47 8.19 7.96 
B 6.010" 4.310" 3.00 4.49 8.09 7.25 
Cc 3.6 X10" 4.6 X10” 3.00 4.50 8.19 7.19 
D 3.010” 9.1108 3.00 4.56 8.33 7.42 
E 2.6 X10" 3.210" 3.00 4.50 8.25 7.47 
F 2.6 X10" 1.110" 3.00 4.50 8.12 7.06 
G 2.1X 10" 4.310"? 3.00 4.49 8.16 7.08 
H 1.810" 4.4107 3.00 4.52 8.26 7.35 
I 7.5 X10" 4.5108 3.00 4.51 8.22 7.14 
x 6.010" 8.3 10!7 3.00 4.50 8.17 6.99 
Kt 6.0 X10" 8.3 X10"? 3.00 4.49 8.05 6.88 








*The yields of Xe'** have not been corrected for neutron capture in Xe, 

{Normalization point. 

{Samples J and K were irradiated at Oak Ridge. The remaining samples were irradi- 
ated in the N.R.X. reactor at Chalk River in different positions. 


mate neutron fluxes at which these samples were irradiated are also shown. 
The large variation in the yield of Xe'* in these samples is the result of varying 
amounts of neutron capture in Xe. (The amount of capture is dependent on 
the length of irradiation as well as the neutron flux.) As can be seen from 
Table II, there are smaller, but still significant, variations in the yields of 
Xe'*4 and some indication of variations in the yields of Xe". These latter two 
effects are not the result of any known or suspected neutron capture process. 
Since some of these samples were irradiated in positions in which a fairly high 
flux of fast neutrons is present, the variations in the yield of Xe" (and Xe! 
if there is a real effect at this mass) may be due to a considerable contribution 
of U8 + n fission in some of the samples. 

The fission product isotope abundance data for the new irradiations are 
given in Table III. The deviations shown in this table are the ‘‘standard 
deviations” of the mass spectrometer measurements. As stated previously, 
irradiation 1 should represent 100% U** + n fission and irradiation 4 100% 
U*8 + nm fission. Irradiations 2 and 3 were carried out under conditions in 
which appreciable contributions of both U** and U8 fission would be expected 
with the larger contribution of U** fission expected in irradiation 3. It can be 
seen from Table III that an increase in the amount of U*§® fission results 
mainly in a decrease in the relative yield of Xe!** accompanied by an increase 
in the relative yields of Xe'*! and Xe'*. In all these later irradiations variations 
in the Xe!* yield due to neutron capture in Xe! are small because of the low 
flux of thermal neutrons. The Xe" yield correction for this process is 0.7% 
for the thermal column irradiation and negligible in other cases. The Xe 
yield data in Table III for the different irradiations are therefore comparable 
and free from neutron capture effects. 

To obtain a comparison between the results of the new irradiations and the 
data given in Table II, the yield data for the four new irradiations have been 
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TABLE III 
MASS SPECTROMETER ABUNDANCES OF FISSION PRODUCT XENON 


Atom per cent 
Irradiation Sample 


Xe! Xels2 Xe Xe36* 

1 Thermal column 1 13.43+.05 20.10+.05 36.78+.04 29.69+.08 
2 13.41+.04 20.03+.05 36.814.05 29.75+.07 
3 13.344.13 19.984.13 36.89+.14 29.79+.20 

Average 13.39 20.04 36.83 29.74 
2 Cd shielded (1) 13.67+.07 20.67+.07 36.214.08  29.45+.06 
(2) 13.714.03 20.36+4.04 36.14+4.07 29.79+.06 

Average 13.69 20.51 36.18 29.62 

3 Los Alamos 

Normal U #1 13.724.08 21.814.09 34.25+.09 30.22+.09 
#2 13.72+.06 21.744+.05 34.20+.05 30.34+.05 

Average 13.72 21.77 34.23 30.28 
4 Depleted UO; 41 15.394.14 22.96+.12 32.714.18  28.94+.09 
#2 16.06+.11 23.154.13 32.16+.09 28.63+.13 

Average 15.72 23.05 32.44 28.79 








*The Xess ‘yield for the thermal column samples should be decreased by 0.7% to correct for 
neutron capture in Xe, This correction should be negligible for the other samples. | 


TABLE IV 
FISSION PRODUCT XENON ABUNDANCE DATA NORMALIZED AT MASS 131 


Relative abundance 





Irradiation Sample ———_——- ——_—_—_—_—_—— 
Xe! Xe! Xe! Xe36 
1 Thermal column 3.00 4.49 8.25 6.62* 
(U*% fission) 
2 Cadmium shielded 3.00 4.49 7.93 6.49 
normal uranium 
3 Los Alamos 3.00 4.76 7.48 6.84 
normal uranium 
4 Los Alamos 3.00 4.40 6.19 5.49 


depleted uranium 
(U8 fission) 





*The yield of Xe'* in the thermal column sample has been corrected for neutron 
capture in Xe}, 


normalized to 3.00 at mass 131 in Table IV. It is obvious that the results for 
the 11 samples in Table II resemble very closely those for the thermal column 
irradiations, indicating that the amount of U*** fission in the former samples is 
not large. 

The data in Table IV should not be interpreted as an indication that the 
absolute fission yield of Xe"! is 3.00% in both U** and U* fission. No highly 
accurate radiochemical yield measurements have been reported for neutron 
fission of U8 in the xenon mass range.* Therefore it is not possible to normalize 
the relative xenon fission yield data obtained for U** fission (depleted uranium 


*Keller, Steinberg, and Glendenin have recently reported a yield of 4.7+0.7% for Te'!® (9). 
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irradiations). Although a comparison of the actual yields of the xenon isotopes 
in the U** and U* fission samples is not possible, several conclusions can be 
drawn from the relative yield data. 

The fission product xenon abundance data for the new irradiations are 
plotted in Fig. 3. The right-hand ordinate scale in this figure gives the fission 
yields for the thermal column irradiations (U** fission). The curve in Fig. 3 


35 







7 
° 


ow 
So 


- ge 
° 
» THERMAL COLUMN IRRADIATION (%) 


XENON ISOTOPE ABUNDANCE (ATOM PERCENT) 


°Q 
a 
w 
> 
20 z 
3 
@! NORMAL U THERMAL COLUMN N-RX- 2 
+2 NORMAL U CADMIUM SHEATHED WRX °° 
O3 NORMAL U FAST REACTOR 
15 @4 DEPLETED U FAST REACTOR 
3-0 


136 


MASS 


Fic. 3. Fission product xenon isotope abundance variations. 


indicates the approximate position of the smooth radiochemical fission yield 
curve for U** fission (4). Since radiochemical studies (4) indicate that the 
smooth yield curve for U** fission parallels that for U** fission in the xenon 
mass range, the curve in Fig. 3 should also represent the position of the U*** 
smooth curve compared to the experimental points for the Los Alamos de- 
pleted uranium irradiation. Therefore, when the relative yields of the xenon 
isotopes are plotted as atom per cent in Fig. 3, it is possible to measure the 
magnitude of fine structure effects from one smooth curve. In this manner the 
variation of fine structure can be studied without the necessity of knowing 
absolute fission yield data for the different irradiations. 

The amount of fine structure (i.e. the percentage increase of the observed 
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vield over the smooth curve value) observed in each of the four different 
irradiations is indicated on Fig. 3. It is obvious that the fine structure at mass 
134 decreases in going from U*** to U8 fission (thermal column irradiation to 
Los Alamos depleted uranium irradiation). At the same time, fine structure 
appears at masses 132 and 131 in U8 fission which was not evident in U?* 
fission. 

In order to interpret the results of the cadmium sheathed and Los Alamos 
normal uranium irradiations, it was first assumed that these samples contained 
contributions of both U?* fission xenon (with the abundances observed for 
the thermal column irradiations) and U8 fission xenon (with the abundances 
observed for irradiation number 4). This assumption, of course, will only be 
completely valid if the relative yields of the xenon isotopes in U?*5 + m fission 
do not change when the neutron energy is altered. With this assumption it is 
possible to calculate the contribution of fission xenon from U*® fission in 
irradiations 2 and 3 as follows: 

Let Y,,55 = observed abundance mass m in U** fission (irradiation 1), 

Y 38 = observed abundance mass m in U8 fission (irradiation 4), 


Ym” = observed abundance mass ™ in irradiation 2 or 3. 


If a is the fractional contribution of xenon from U8 fission in irradiation 2 
or 3, then 
Ym? = (1—a) Yp28+a Y,,78. 


The values of a found by substituting the abundance values from Table III 
in this equation are given in Table V. The contributions of xenon from U8 


TABLE V 


CALCULATED CONTRIBUTION OF XENON FROM U8 FISSION TO THE 
OBSERVED ABUNDANCES IN IRRADIATIONS 2 AND 3 





Fractional contribution of xenon from U*® fission 
Calculated from mass_=——— 


Irradiation 2 Irradiation 3 
131 13 14 
132 16 .58 
134 15 59 
136 13 —.57 
Average .14 


fission calculated from the different xenon isotope abundances for irradiation 2 
agree very well. The calculations indicate that for the cadmium sheathed 
irradiation 14% of the fission product xenon resulted from U** fission. However, 
the calculated U* fission contributions for irradiation 3 from the different 
xenon isotope abundances do not agree with each other. This disagreement is 
almost certainly due to the fact that the U* fission yields at high neutron 
energies are not the same as those observed with thermal neutrons. 

The samples for irradiation number 2 were normal uranium sheathed in 
cadmium and were irradiated in the N.R.X. reactor at Chalk River. The 
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cadmium sheath on these samples would absorb all thermal neutrons but the 
absorption of neutrons with energies greater than ~ 1 ev. would be small. 
Since there is a considerable epithermal neutron flux at the position in the 
N.R.X. reactor at which these samples were irradiated in addition to the flux 
of fast ‘‘fission spectrum’’ neutrons, the majority of the U™* fission in these 
samples would be produced by epithermal neutrons and not by the fast 
neutrons. Therefore the effective neutron energy for the U*™® fission contri- 
bution in irradiation 2 was probably not very high. For this reason the U** 
fission contribution in irradiation 2 would be expected to produce fission xenon 
with essentially the same isotope abundance as observed in thermal neutron 
fission of U**. In this case, therefore, the assumption made in calculating the 
U*8 fission contribution is valid and the good agreement shown in Table V, 
column 2, would be expected. 

The samples used in irradiation 3 were normal uranium and were irradiated 
in the fast reactor at Los Alamos. In this reactor the neutron energy distri- 
bution is essentially the fission neutron spectrum and therefore the U™® fission 
in these samples will be produced by high energy neutrons. Therefore the 
assumption made in calculating the U*** fission contribution for this irradiation 
(Table V, column 3) will not be valid unless the U* fission yields do not 
change with relatively large changes in neutron energy. Since it is known (12) 
that the fission yields do vary when the neutron energy is changed by a large 
amount, the disagreement in Table V, column 3, is not surprising. 

There is good agreement between the calculated U** fission contributions 
for irradiation 3 at masses 132 and 134. Also the difference between the U*® 
and U8 yields is greatest at mass 134, and thus less likely to be greatly altered 
by neutron energy changes. For these reasons, it is reasonable to assume that 
the contribution of fission xenon from U** fission is about 59% for irradiation 3 
and that the discrepancies at mass 131 and 136 in Table V, column 3, are the 
result of the variation with neutron energy of the xenon yields from U** fission. 

It should be noted that the 14% and 59% contributions of fission xenon 
from U** fission in irradiations 2 and 3 respectively do not necessarily indicate 
that U8 fission accounted for 14% and 59% of the fission in these samples. 
The latter is only true if the absolute fission yield of the stable xenon isotopes 
is the same in both U®* and U** fission. Actually radiochemical yield studies (4) 
indicate that the smooth yield curves for U?* and U*® fission are quite close 
together in the mass range 131-136. Therefore the total yield of xenon in each 
type of fission is probably nearly the same. On this basis it can be said that 
for the cadmium sheathed irradiation approximately 15% of the fission was 
U8 + n fission. The corresponding figure for the normal uranium irradiated 
at Los Alamos is about 60%. Furthermore it can be seen from Fig. 3 that, if 
the total yield of xenon is the same for U™® and U*** fission, the fission yield of 
Xe!*! in U*8 fission (depleted uranium irradiation 4) is approximately 3.50% 
compared to 3.00% for the same isotope in U** fission. 

The ‘‘neutron boil-off’” mechanism suggested by Glendenin and elaborated 
by Pappas to explain the fine structure observed in U** + n fission would pre- 
dict a higher yield at mass 131 and 132 in U™® fission than in U** fission. This 
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is in agreement with the results given in Table III. A quantitative comparison 
between the predicted and experimental yields is given in Table VI where the 


TABLE VI 


PREDICTED AND OBSERVED FINE STRUCTURE EFFECTS 
IN XENON FISSION YIELDS 











235 238 
Mass a a 
Predicted Observed Predicted Observed 
131 3.00* 3.00* 3.50* 3.50* 
132 4.65 4.49 5.10 5.04 
134 7.74 8.21 7.31 7.01 
136 6.50 6.48 6.63 6.24 


*Normalization point. 


predicted yields have been calculated following the method adopted by 
Pappas (11). In this table the U** results have been arbitrarily normalized at 
3.50% for Xe'!. It can be seen that the quantitative agreement is quite good. 
However, it should be noted that the predicted fine structure is too low for 
U* fission and too high for U** fission. Recently (5) it has been found that the 
discrepancy between the predicted and experimental yields is much larger 
in the case of U** fission. This indicates that while the mechanism developed 
by Pappas offers a reasonably accurate explanation of the observed fine 
structure in the xenon range in U** and U*%8 fission (Table VI), it must be 
further modified before it will completely explain all available experimental 
data. 
(ii) Krypton Isotope Fission Yields 

Table VII contains the krypton isotope abundance data for the irradiations 
listed in Table I. In normalizing the krypton results a value of 5.05 was 


TABLE VII 
MASS SPECTROMETER ABUNDANCES OF FISSION PRODUCT KRYPTON 





Atom per cent 





Irradiation Sample 
Kr83 Kr Kr8* Kr'6 

1 Thermal column 1 14.16+.06 26.32+.13 7.50+.05 52.02+.13 
2 14.07+.06 26.10+.14 7.55+.06 52.28+.07 
3 14.074.07 25.388+4.14 7.724.15 52.83+.25 

Average 14.10 25.93 7.59 52.38 
2 Cd sheathed (1) 14.22+.01 27.09+.05 7.31+.06 51.38+.09 
(2) 14.044.13 26.774.07 7.67+.09 51.52+.16 

Average 14.13 26.93 7.49 51.45 

3 Los Alamos 

Normal uranium #1 15.08+.30 28.52+.11 6.934.20 49.47+.28 
#2 14.85+.38 27.62+.42 6564.24 50.97+.42 

Average 14.97 28.07 6.74 50.22 
4 Depleted UO; #2 14.50+.60 30.50+1.40 5.50+.50 49.50+1.10 





*Kr® corrected assuming T3 = 10.27 years. 








FISSIOM YIELD (%) 


WANLESS AND THODE: FISSION YIELDS 





20 


Fic. 4. Kr® decay scheme (7). 


@1! NORMAL U THERMAL COLUMN N-R:X- 
+ 2NORMAL U CADMIUM SHEATHED N-R:X: 
O3 NORMAL U FAST REACTOR 


@4 DEPLETED U FAST REACTOR 


84 85 
MASS 


Fic. 5. Krypton fission yields. 





86 


51 











552 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


emploved for the ratio of fission product xenon to fission product krypton in 
U*5 fission (1). Since there is a discrepancy between reported values of this 
ratio (1, 13), a precise determination is at present under way in this laboratory. 
The normalization employed however will not alter the conclusions drawn here. 

The krypton isotope abundances for irradiation 1 (U** + n fission) have 
been normalized at 0.62% for mass 83. The krypton isotope at mass 85 has 
two isomeric states (Fig. 4), a long-lived one with a half-life of 10.27 years (14) 
and a short-lived one with a 4.4-hour half-life (7). Since the shorter-lived 
isomer will have completely decayed before the fission gases were extracted 
from the uranium samples, the mass spectrometer results will indicate only 
the yield of the 10.27-year isomer. All Kr® yields have been corrected for the 
decay of the long-lived isomer up to the time of analysis. In order to compare 
the fission yields at this mass, the Kr*® yield for irradiation 1 has been arbi- 
trarily placed on the smooth curve (see Fig. 5) and the yields for the other 
samples have been increased proportionately. This corresponds to a branching 
ratio for the decay of Kr® of 0.29. The value of this branching ratio measured 
directly from 6 decay data is 0.28 (3). The assumption that the mass 85 yield 
falls on the smooth curve for U** fission (sample 1) would therefore appear to 
be justified. The 85 mass chain yields for the other samples have been adjusted 
assuming that the branching ratio remains constant as the fissioning nucleus 
is changed from U** to U*8, This latter procedure is justified since, on the 
basis of Glendenin’s charge distribution theory, the primary yield of Kr** and 
Kr®" is negligible in both cases (6). 

There is a decided shift of the light mass peak of the fission yield curve 
between U** and U*® fission (4). In Fig. 5 the depleted UO; yields (irradi- 
ation 4) have been displaced at mass 86 by 0.86% (see Table VIII) to corre- 

TABLE VIII 
KRYPTON FISSION YIELDS 


% Fission yield 


Irradiation Sample 
Kr® Kr* Kr®* Kr86 
1 Thermal column 0.62 1.14 1.48 2.30 
2 Cd sheathed 0.60 1.14 1.41 2.18 
3 Los Alamos 0.53 0.99 1.06 1.78 
normal uranium 
4 Los Alamos 0.42 0.89 0.68 1.44 


depleted uranium 


*The mass 85 chain has been corrected for the decay of the 4.4-hour 
isomer of Kr® by arbitrarily placing the mass 85 yield on the smooth mass— 
yield curve for sample 1 (U*® + n fission). This corresponds to a branching 
ratio for the decay of Kr® of 0.29. The value of this branching ratio meas- 
ured directly from B decay data is 0.28 (3). The assumption that the mass 
85 yield falls on the smooth curve for U** fission (sample 1) would therefore 
appear to be justified. The 85 mass chain yields for the other samples have 
been adjusted assuming that the branching ratio remains constant as the 
fissioning nucleus is changed from U** to U*8, This latter procedure 1s 
justified since the primary yield of mass 85 is negligible in both cases. 


spond to the above-mentioned shift. The Los Alamos normal uranium results 
(irradiation 3) have been adjusted at mass 86 in accordance with a 60% 
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U*8 + n fission in this sample (as indicated from the xenon data). Finally, 
the cadmium sheathed results (irradiation 2) have been incorporated in Fig. 5 
by adjusting the mass 86 yield to correspond to 15% U*** + n fission as de- 
termined from the xenon results. It is now possible to draw four roughly 
parallel smooth curves from which the percentage variations may be calcu- 
lated. Again if fine structure is defined as the percentage variation from the 
respective smooth curves, these results indicate large fine structure shifts at 
masses 84 and 85. 

It will be noted that at mass 84 the fission yields are all higher than the 
corresponding smooth curve values while at mass 85 the opposite is found to 
be the case. The fine structure shift in going from U?* + n fission to U8 + n 
fission amounts to 388—20 = 18% at mass 84 and 29% at mass 85 (see Fig. 5). 

The fine structures determined for the Los Alamos normal uranium and the 
Chalk River cadmium sheathed irradiations, 3 and 2 respectively, fall in the 
intermediate range and indicate 10/18 = 56% and 4/18 = 22% U** +n 
fission respectively when calculated from the mass 84 fission yields. However, 
when calculated from the yield values at mass 85, the corresponding U** + n 
fission contributions indicated are 52% and 0%. The U** + n contribution 
indicated for the Los Alamos normal uranium sample is in good agreement 
with the value found in the xenon range (60%). In the case of the cadmium 
sheathed sample the 22% U** + n fission contribution indicated at mass 84 
agrees with that found in the xenon range at mass 134, but does not agree at 
mass 85. This latter discrepancy may be due to the difficulty experienced in 
measuring the small peak recorded at mass 85 (see Fig. 1). 

The occurrence of fine structure in the krypton mass range has now been 
definitely established for the first time. This fine structure changes in going 
from U* + » to U8 + n fission. Specifically, the yield at mass 84 increases 
while the yield at mass 85 decreases for U*** fission compared to the values 
obtained for U** thermal neutron fission. The fission yields predicted by the 
method of Pappas have been calculated and no agreement has been found with 
the experimentally determined yields. 

Since the observed fine structure does not appear to be the result of a 
Glendenin type ‘‘neutron boil-off’” from primary fission fragments with 51, 
53, 55, and 57 neutrons, some other explanation must be found for the observed 
fine structure. It is suggested that a delayed neutron emitter in the 85 mass 
chain produces the observed fine structure at masses 84 and 85. The tentative 
assignment of such an activity to As® (the actual neutron emission occurring 
after the 6-decay of As*) is possible on the basis of nuclear systematics. Since 
the cumulative yield of As* would be expected to be higher for U*** fission than 
for U** fission on the basis of Glendenin’s charge distribution theory (6, 11), 
the fine structure resulting from delayed neutron emission by As* would be 
more pronounced in U* fission. This is in qualitative agreement* with the 
experimental results (Fig. 5). 

*It has been shown that a reasonable quantitative explanation of the observed fine structure in 
the krypton region can be obtained not only for U* and U™ fission but also for U** and Pu®® 


jission and for the spontaneous fission of U™* by assigning the 0.48 sec. delayed neutron activity 
to As® (W. H. Fleming, Ph.D. thesis, McMaster University, 1954). 
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THE PROTON COMPONENT OF COSMIC RAYS AT 
AIRCRAFT ALTITUDES! 


By K. W. OGILviE? 


ABSTRACT 
Experiments with a heavy particle selector installed in an aircraft show 
that the proton intensity in the energy range 230 to 480 Mev. has a vari- 
ation with height which is exponential, with absorption length 128+10 gm./ 
cm.? The latitude effect between A = 40°N. and A = 56°N. is in general agree- 
ment with that found for fast (up to 5 Mev.) neutrons. Some information has 
been obtained on the variation of the average specific yield for primary energies 
below 5 Bev. 
INTRODUCTION 
The altitude and latitude dependence of the nucleonic component of the 
cosmic radiation have been studied (13) using neutron counters carried in 
aircraft. The present experiments supplement this work in the region between 
geomagnetic latitudes 56° and 40°, and between ground level and 25,000 ft., 
by using a ‘heavy particle detector’ that records the passage of protons of 
energy between 230 and 480 Mev. 
APPARATUS 
The apparatus, Fig. 1, consists of a Cerenkov counter of large area, 16 in. 


X 16 in., containing water, placed in a Geiger counter telescope. Passage of 
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Fic. 1. Block diagram of apparatus. Anticoincidences between the outputs of the top two 
coincidence circuits denote heavy particles; anticoincidence circuit not shown. 


a particle through the telescope, without a corresponding count in the Cerenkov 
counter, indicates that the particle was travelling at a velocity below the 
Cerenkov threshold. This quantity is V. = c/m, where c is the velocity of 
light in free space and m is the refractive index of the substance used as the 
Cerenkov radiator. Deurden and Hyams (2) have shown that absorbers can 
be placed in the Geiger telescope in such a way that all the particles with mass 
less than a certain figure (u mesons) that penetrate them are travelling faster 
than V,, and that in this condition the apparatus counts particles with greater 


1Manuscript received June 17, 1956. 

Contribution from the Division of Pure Physics, National Research Council, Ottawa, Canada. 
Issued as N.R.C. No. 8707. 

2National Research Laboratories Postdoctorate Fellow. 
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mass (protons) in an appropriate energy range, by means of recording anti- 
coincidences between the Geiger telescope and the Cerenkov counter. A rela- 
tivistic particle must travel through the tank for a distance less than 20% of 
its depth to be recorded. A more complete description of the apparatus will 
be contained in a paper to be published shortly. With the present apparatus 
with an absorber of 3.7 cm. Pb, sensitivity to protons is obtained in the energy 
range 230 to 480 Mev., in which region the intensity at sea level is 3.03 X10~° 
particles/cm.?/steradian/sec. (6). 
RESULTS 


The Altitude Dependence of the Intensity of Protons 

The apparatus, mounted in two racks, was installed in a ‘North Star’ air- 
craft of the Central Experimental and Proving Establishment, R.C.A.F., 
Rockcliffe. The aircraft was then flown at 5000-ft. levels up to and including 
25,000 ft. for periods ranging from three hours at 5000 ft. to one hour at 
25,000 ft. Altitude was determined by the pilot’s altimeter, which was set at 
29.92 in. on the ground before each flight, and checked by a similar instrument 
mounted close to the apparatus. These altimeters are related to the N.A.C.A. 
standard atmosphere, which was used in the reduction from altitude to pressure 
in gm./cm.? (7), and fluctuations in altitude were less than 100 ft. at 25,000 ft. 
The flights followed courses closely parallel to geographic latitude 45° 30’, 
which at this longitude is roughly parallel to geomagnetic latitude 57° and, 
in any case, is close to the knee of the latitude effect. The courses followed at 
the different levels during flight 1 are shown in Fig. 2. The results obtained 
for the ‘heavy particle’ rate are set out in Table I, and those for the ‘hard 
component’ in Table II. 





Fic. 2. Courses for constant altitude run 1. 


Fig. 3 shows these data plotted on a semilog scale, together with the best-fit 
straight line for the data of Table I. This gives the value 128+12 gm./cm.? 
for the absorption length of the particles observed, which are to a good approxi- 
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TABLE I 
VARIATION OF ANTICOINCIDENCE COUNTING RATE WITH ALTITUDE 





Flight and date 
Altitude, ——————_——————————— SOU Weighted 














ft. 1 2 3 4 mean 
17/3/55 24/3/55 25/3/55 28/3/55 

25,000 8832 +150 7500 +200 8160+150 
20,000 3800 + 100 3400+ 100 3600 +70 

15,000 2160+100 1500 +80 2010+60 1890 +60 

10,000 576+20 546+20 560+20 

5,000 272+20 175+12 194+15 
Sea level 516+0.5 
N.B.—AlI1I errors statistical. 

TABLE II 
VARIATION OF TRIPLE COINCIDENCE COUNTING RATE WITH ALTITUDE 
Flight and date 
Altitude, —— Mean 
ft; 1 2 3 4 
17/3/55 24/3/55 25/3/55 28/3/55 

25,000 136800 137500 137200 
20,000 90000 92700 91300 

15,000 59400 54300 60500 58070 

10,000 38900 38400 38650 

5,000 26600 26500 26500 
Sea level 19300 





Examination of the records of the neutron recording pile at Ottawa for the period 9 to 31 March, 
1955, disclose an intensity increase of 3$% on the 13 which had decayed to 14% by the 17 March. 
Small fluctuations of the order of 1% persisted through the period of the experiments. 
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Fic. 3. A. Altitude variation of proton rate. B. Altitude variation of triple coincidence 
rate. 
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mation only protons. The good straight line obtained for the proton rate 
variation is inconsistent with a significant proportion of the anticoincidence 
counts being due to slow yu-mesons, taking into account the much smaller 
increase of the u-meson rate with altitude. Some proton counts could be due 
to deuterons, since the altitude variation of the deuteron intensity is likely to 
be similar to that of protons, but the proportion is thought to be less than 10%. 
The lower mass limit of particles to which the apparatus is sensitive is approxi- 
mately 550 electron masses, so that m-mesons cannot be detected. 

As may be seen from Fig. 1, an accidental triple coincidence in the telescope 
will be recorded as a proton, if unaccompanied by a cancelling pulse from the 
Cerenkov counter. This point has been investigated in connection with other 
experiments with an identical apparatus at sea level. The total rate of triple 
coincidences is given by the expression 


Ne = N+3NiN2N3T?+2Mi2N3T+2Ni3NoT +2Na3MiT, 


where JN is the number of real triple coincidences, N, the counting rate in the 
rth counter, and N,, the actual coincidence rate between counters 7 and s. 
With the present value of the resolving time, T = 2 usec., the contribution of 
the second term may be neglected; the last two terms involve the traversal 
of the tank by relativistic particles, and thus need not be considered. The 
third term then remains, and must however represent a large overestimate of 
the actual number of spurious protons. This is because many of the particles 
adding to the term 2N,.N3T, by their traversals of counter trays 2 and 3, also 
go through enough of the tank to be detected and give rise to a cancelling 
pulse. The term 2Ni2N3T was measured during the experiments, using the 
circuit shown dotted in Fig. 1, and remained throughout at about 10% of the 
total; no correction has been applied for this effect. 

In order to convert the counting rates to absolute intensities it is necessary 
to evaluate the solid angle of the apparatus using the triple rate at sea level 
and the known hard component intensity. The values given in Table III have 


TABLE III 
PRESSURE VARIATIONS OF THE ABSOLUTE INTENSITIES (IN PARTICLES/CM.?/SEC./STERADIAN) 





_ Proton Hard component Hard component 
Pressure intensity intensity (7.6 cm.) intensity (10.4 cm.) (9) 


383 508+10 x 1075 6.12 X 107 5.2 X 107? 
475 24445 4.07 3. 

583 117.544 2.59 

710 34.84+1.2 LIZ 

860 10.2+0.8 1.18 

1030 3.2+0.3 0.86 (assumed) 


been obtained in this way, and have also been corrected for the loss of particles 
from the lower end of the energy interval due to the lengthening of the path 
through the absorber by scattering. This correction, amounting to 19+5%, 
is evaluated in the following way. The lateral distribution about the pro- 
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jection on the lower counter tray of their point of entry to the absorber may 
be calculated for protons of various energies. 

Assuming their paths to be initially vertical, there will be, for each energy, 
numbers of protons unable to penetrate the absorber, namely those for which 
the deflection due to scattering is greater than x, given by 


x? = (0.907R)?—d?, 


where d is the absorber thickness, here 3.7 cm. 

The factor 0.907 = T/R is given by Koenig (3), where J is the straight 
track in an absorber corresponding to an actual path of length R followed by 
the particle. At 270 Mev., for example, 40% of the protons will not penetrate 
to the lower counter, while at 320 Mev. only 63% will fail to be recorded. In 
order to work out the correction the shape of the spectrum must be known, and 
it is by taking the spectrum of Mylroi and Wilson that we arrive at the value 
given above. It must be emphasized that the corrections affect the absolute 
values of the proton intensity only, and have an effect on the slope of the 
proton curve that is less than the experimental error. The values shown in 
column four of Table III were obtained by Neher in airplane flights at latitude 
53° N. geomagnetic. 

Discussion 


Simpson (12) finds for the exponential absorption length of the neutron 
producing radiation in the atmosphere at geomagnetic latitude 53° the value 


L = 1538+2 gm./cm? 


The detector used by Simpson was a BF; proportional counter surrounded 
by paraffin, with a cadmium shield, and responded to neutrons with energies 
between 0.4 ev. and about 5 Mev. The atmospheric absorption of high energy 
‘N rays’ is exponential with an absorption length of about 125 gm. cm.? (10). 

Our result, corresponding to a mean energy of about 400 Mev., because of 
the loss of low energy particles by scattering, agrees well with the latter value. 
Rossi (11), using the production spectrum of protons at 50 gm./cm.? measured 
by Camerini ef al. (1), finds the following expression for the intensity of 
protons of energy E at depth x gm./cm.? 


age” En—E 
K(E) (50+En)(50+£)’ 


where E is defined by the relation 
R(Em) = R(E)+<x, R(E) being the range of a proton of energy E. 
Thus 


J’ (E,x) = 


- ==) . LE, i} End 50+Em | 
NL ~ JP(E, x2)LEm,—E 50+Em, 4° 
Working out the bracketed expression, Y, with E = 400 Mev., and x, 
= 1030 gm./cm.*, we arrive at the factors given in Table IV. The quantities 
in column three are plotted in Fig. 4 and yield a value of L = 122412 gm./cm.? 
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Fic. 4. A. Altitude variation of proton intensity. B. Altitude variation of triple coincidence 
rate. 


TABLE IV 
CORRECTED PROTON RATE 


x ¥ JP(E, x) 
gm./cm. 

383 1.26 10280 
475 1.18 4200 
583 1.12 2140 
710 1.07 600 
860 1.03 200.5 
1030 1.00 51.6 


The Latitude Effect 

A flight was undertaken by the same aircraft, carrying the same apparatus, 
from Ottawa along a line joining Ottawa to Nassau, and terminating at Miami. 
The results set out below in Table V show that this course follows longitude 
76° W. which corresponds closely to a constant geomagnetic longitude. The 
altitude chosen was 20,000 ft. (= 475 gm./cm.?), a compromise between 
counting rate and convenience, as the aircraft employed was unpressurized. 
The reduction from geographic to geomagnetic latitude was made using the 
curves of McNish (5). The statistical errors are seen in Table V to be approxi- 
mately +2%, the counting intervals being approximately 25 min. in duration. 
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TABLE V 
RESULTS OF LATITUDE VARIATION EXPERIMENT 





Mean geographic position Mean Counting rate Triple 
geomagnetic of protons coincidence 
Latitude Longitude latitude per hour rate per hour 
Going south 
43° 50’ N. 76° 4’ W. 55° 30’ 4180+80 90,500 +400 
42° 13’ 76° 2’ 53° 53’ 4305 +80 91,600+400 
40° 15’ 76° 9’ 51° 55’ 3980 +80 90,400 +400 
38° 18’ 76° 29’ 49° 58’ 3900 +80 90,700 +400 
36° 7’ 76° 40’ 47° 47’ 3690 +80 88,300 +400 
34° 2’ 76° 50’ 45° 42’ 3900 +80 89,200 +400 
32° 12’ 76° 44’ 43° 52’ 3520 +80 86,500 +400 
30° 35’ 76° 44’ 42° 15’ 3710+80 87,000 +400 
29° 6’ 76° 50’ 40° 46’ 3140+80 83,000 +400 
27° 40’ 76° 43’ 39° 20’ 2945 +80 83,000 +400 
Going north 
33° 42’ 77° 0’ 45° 22’ 3790+80 87,000 +400 
36° 5’ 76° 45’ 47° 45’ 3940+80 88,400 +400 
37° 20’ 76° 30’ 49° 0’ 4000 +80 88,400 +400 
40° 9’ 76° 28’ 51° 49’ 4160+80 _ 
42° 5’ 76° 20’ 53° 45’ 4340 +80 89,500 +400 
44° 14’ 75° 57’ 55° 54’ 4470+80 89,500 +400 





The aircraft position was obtained at frequent intervals, and the mean counting 
rate for each counting interval plotted at the position corresponding to the 
center of that interval. 

We must now consider the possible influence of the longitude effect. In the 
region of longitude 95°-110° W., and at geomagnetic latitude 41°, [d7/d¢], 
= —1% per degree (12), at atmospheric depth 448 gm./cm.*, for neutrons. 
Reference to the table shows that error due to longitude variation during the 
flight is less than the statistical errors. We can compare the slope of the curves 
and the magnitude of the effects for neutrons and protons, but not the in- 
tensities, by assuming that the variation of the longitude effect with latitude 
is the same in the region of longitude 76° W. as in the region of 116° W. In 
Fig. 5 our results are plotted together with those of Simpson (12), curve B, 
for fast neutrons at 448 gm./cm.? (23,000 ft.) and longitude 116° W. Corre- 
spondence is close, but our results suggest a knee at higher latitudes than the 
neutron results. The accuracy is hardly high enough to settle this, however. 
Table V also contains details of the triple coincidence rates recorded at the 
same time as the proton counts. These are plotted in Fig. 6. A knee can be 
observed in this curve at \ = 53°, and the observed intensity difference found 
between A = 54° and A = 39° is 8.342%. The total absorber in the telescope 
was 3.7 cm. of lead, 20 cm. of water, and the counter walls. 


The Average Specific Yield 

Below about 200 gm./cm.? in the atmosphere, protons in our present energy 
range observed at depth x gm./cm.? are produced close to this depth. Thus a 
measurement of the vertical proton flux amounts to a measurement of the 
proton production rate at the same depth. The total production rate at depth x 
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Fic. 5. A. Latitude variation of proton intensity. B. Neutron curve due to Simpson. 
Fic. 6. Latitude variation of triple coincidence rate. 


and latitude \, due to primaries arriving into unit solid angle in a vertical 
direction, is given, (14), by the expression 


Rie) = — J” SAN, x) fa )dFr/AN] AN, 


where dF,/dN is the total differential flux of primary particles per unit solid 
angle having momentum to charge ratio equal to N, fa(N) is the fraction of 
primaries having atomic weight A, and S,(N, x) is the specific yield. This 
latter quantity is, in our case, the production rate of protons in our energy 
range, due to a unit vertical flux of primary particles of atomic weight A and 
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momentum-:to charge ratio NV. The average specific yield is then given by the 
relation 


S(N,x) = » Sa(N, x) fa(N) 


so that 
S(N, x) = dR, (A, x)/dF7(N). 


The quantity R, may be equated to the counting rate of our apparatus, 
corrected to absolute intensity. The latter procedure entails errors of less than 
20%, and the Gross transformation need not be applied to our results, contrary 
to the procedure adopted by Treiman (14), since a telescope is used here. In 
the case of exponential variation of the counting rate with altitude, the Gross 
transformation reduces to multiplication by a factor [1+x/L]/2z, where L is 
a function of }. Between 55° and 20° latitude, Z varies from 125 to 135 
gm./cm.? (4). We can obtain the vertical cutoff values of momentum to charge 
ratio, (A), for the values of \ at which our observations were made, from the 
results of Neher (8), and the corresponding quantities F (JV) from the spectrum 
of Winckler (15). Following Treiman, we plot the vertical proton production 
rate at 475 gm./cm.? against the vertical primary flux, Fig. 7. On differentiating 
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Fic. 7. Vertical intensity of protons plot- Fic. 8. Yield of particles in our energy 
ted against vertical primary flux. range as a function of P/Z. 


the smooth curve shown, we obtain the relation between the average specific 
yield and the momentum to charge ratio, Fig. 8. It will be seen that observa- 
tions were not carried out sufficiently close to the equator to show any levelling 
off of the yield curve, which is still rising steeply at a primary particle (proton) 
momentum of 5 Bev./c. A primary proton can just penetrate 475 gm./cm.? of 











564 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


air when its momentum is 2.1 Bev./c, so that a primary proton can be detected 
by our apparatus when it has a momentum of 2.8 Bev./c. Neglecting the 
possibility of a primary particle arriving at our level by travelling a large 
portion of the range as a neutron, the yield at P/Z = 2.8 Bev./c should be 


Y(P/Z = 2.8) = e475 /L6 


where L; is the interaction length of such a proton. This follows from the fact 
that any interaction suffered will cause a loss of energy and so eliminate the 
particle from our counting rate. The relation leads to a value L; = 71 gm./cm.?, 
this being regarded more as a check upon the accuracy of the absolute intensi- 
ties deduced from the counting rates, rather than as a determination of the 
interaction length of high energy protons. 


Conclusions 

From the foregoing we may conclude that the altitude dependence of 
400 Mev. protons is exponential with an absorption length of 122+12 gm./cm.’, 
taking into account the modification of the spectrum with height. The latitude 
effect for such protons shows a similar variation to that of fast neutrons 
between \ = 56° and A = 40° geomagnetic, and the amount of the effect is 
28.4%, but there is some evidence that the knee of the latitude effect is above 
A = 54°. 
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FIELD EQUATIONS FROM PARTICLE EQUATIONS! 


By W. J. ARCHIBALD 


ABSTRACT 
A simple point of view is developed which shows how various types of field 
equations may be obtained from the relativistic Hamiltonian of a particle. 

The energy in the various fields (electromagnetic, meson, etc.) is concen- 
trated in small regions of space and as these concentrated bundles of energy 
are moving rapidly, they appear as particles to an observer. Sometimes these 
bundles of energy have zero inertial mass (photons) and sometimes they do 
not (mesons). This suggests that there must be some simple connection be- 
tween field equations and the equations of motion of a particle. Of course, 
the quantization of the field equations explains satisfactorily the existence of 
the ‘‘photons in the field’. What is attempted here is to show that, without 
resorting to field quantization, field equations and particle equations are 
intimately related. 


ELECTROMAGNETIC FIELD 


As a simple example consider Maxwell’s equations for the electromagnetic 
field in empty space. It is sufficient to take 


(1] VXE=—-H/, VXH=E/c, 

since V - E = V - H = 0 follow from these if they are imposed as conditions 
at ¢ = 0. From [1], the equation for F = E+7H is 

[2] F/c = —i¥ XF. 

On multiplying by zh this may be written 

[3] OF = icpxF, 


where § — ihd/dt and p— —ih¥V. Equation [3] can be written in matrix 
form as 














| F, | 0 a Py | F, 
[4] © | Fy| = ic| pz 0 ~ bs || Fy 
| F. | | —Py pz 0 || F: 
If only the operator part of [4] is retained, it becomes 
[5] H = os p ’ 
where 
0 0 0 0 0 0-i O 
[6] tz =|0 O-+), Ty = oOo O -O1, ga=lsé €& © 
0 « O —1 0 0 oO oO @ 





Equation [5] looks like the Dirac operator for a particle of zero rest mass, and 


\Manuscript received May 20, 1956. 
Contribution from the Department of Physics, Dalhousie University, Halifax, Nova Scotia. 
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the well-behaved solutions of [4] would certainly be consistent with physically 
significant solutions of [1]. 

If [5] is accepted as the Hamiltonian operator for a particle, then certain 
properties of this particle follow immediately: e.g. the operator M,+Ar,, 
where M, = xpy—ypz, commutes with § so that it is a constant of the motion. 
Thus hr, is the operator associated with the z component of the particle’s spin, 
and this operator has the eigenvalues h, 0, and —f. In confirmation of this 
picture it is easy to show that *«X* = it, which must be true if fr is to be a 
spin operator. Furthermore, | y |? = | Fz |?+ | F, ?+| F. |? = #2?+H?. Thus 
the probability of finding the particle in a volume of given size is greatest in 
those regions where the classical field energy density is greatest. In summary, 
this particle has unit spin and is most likely to be found where the fields are 





strongest. 


MESON FIELD 


For the second example the process is reversed: we start from the equation 
for a particle and arrive at the equations for a field. The relativistic energy 
equation for a free particle of rest mass mp and energy E may be written as 


(imyc)? = ((E/c)P?+p/2+p/+p2, 


or 
[7] (tmoc)? = pe+pr+p2+p3 = P*, 

where obvious abbreviations have been introduced. The four p’s are the 
components of a 4-vector of constant length. To set up the quantum mechan- 
ical equation for a particle Dirac used a symbolic 4-vector (with components 
Qo, @, @2, a3) having the property that (ap)? = P*. (a-p is the four dimen- 
sional scalar product). He then quantized the equation imoc = a-p. This 
particle turned out to be the electron. Let us, without giving any a priori 
reason for doing so, seek a symbolic 4-vector such that 


[8] (a-p)* = P*(a-p). 

Needless to say, the Dirac a’s satisfy this equation, but there may be other 
a’s which do so also. If [8] is to hold the a’s must have the following properties: 
a? =a, (u = 0, 1, 2, 3), 

[9] Ayaytayaa,taa’” =a, (u ¥ »), 
(an pOy + ara yar) + (ayanar+ayanay) + (a,ara,+ta,a,a,) = 0 (A, u, v different). 
While these equations are satisfied by the Dirac operators, they can also be 
met by the following scheme: 
a? =a, (u = 0,1, 2,3), 
[10] avataa,” =a, (ux), 
ayaa, =O (u ¥»v), 


MAyaertaa,on =O (A, uw, v different) . 
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The algebra defined by [10] has been studied by Kemmer (1) and he has shown 
that a’s which are 10 by 10 matrices can be found which obey these equations. 

It is not impossibly difficult to find four such a’s. Since a,’ = a, is the 
simplest equation which each a obeys, each ahas the eigenvalues 1, 0, and —1, 
and if a representation is chosen in which ap is diagonal it seems reasonable to 
have the numbers 1, 1, 1, 0,0, 0,0, —1, —1, —1 along the diagonal. With this 
start it is possible, but tedious, to find three other a’s, and only the results 
will be recorded. Once the a’s are obtained it is easy to show that they have 
the required properties. 

Following Dirac, by analogy, we obtain 


[11] ImMoC = aoPotaipitarp2tasp;, 
and combining the four matrices on the right of [11] 
po A 0 
[12] imypc=|A* O BY, 
0 B —Po 


where (a) each fp is a diagonal 3 by 3 matrix, 
(6) the central zero is a 4 by 4 null matrix, 
(c) the other zeros are 3 by 3 null matrices, 
(d) A and B are 4 by 3 matrices and A* and B* their complex 
conjugates, with 


1 —f, Rs 0 <i. 1 Py Dp: 0 —pz 
({13] A = a Pz 0 -~P: —hy|) B= Fe —Pz 0 Pz —Py!- 
2 0 Pz Py ~f. v2 0 ~Pfy ~“. Pel 


Equation [12] becomes an operator equation by replacing po by —(h/c)d/dt, 
pz by —th d/dx, etc. In this treatment it is not possible to have po multiplied 
by a unit matrix (as is the case for a Dirac particle) since a has no reciprocal. 
The matrices A and B have the following simple properties: 
AA* = BB* = }p’ 1(3), 
[14] AB*B = }p°A, 
BA*A = 3p°B, 
where 1(3) is the unit 3 by 3 matrix and p? = p,?+p,?+).2’. 
It is important to remember that the matrix in [12] is put forward in the 
present instance only because it has the mathematical property that when it is 
cubed it produces itself multiplied by P?. This is certainly not a sufficient 


reason for studying it, but a brief examination discloses interesting properties 
which may have value. Let [12] act on a 10 element column matrix 


v1 
X |» 
ve 


where y, ¥2, and x are three and four element column matrices respectively. 
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Then 
(imoc —po)ti = Ax, 
[15] 1MoCx = A*y,+Br*yo, 
(tmoc + po) v2 = Bx. 
It is difficult to sense the properties of these equations in their present form. 
However it is possible by elimination of the y’s to obtain a second-order 
differential equation for x, which is 
[16] —(mrc?+po)x = (A*A+B*B)x. 


To display all the features of this equation it is necessary to write it out in 
full as 








Ju] |b2+h? pbs xb 0 | x 
17 eB e 2 2 | X2| a Pybz prtp? Prxby : 0 | X2 | 
[ 7] (m Cc + po’) x | —Pprp: Prby pyt+pr 0 X3 
x4 | 0 0 0 Pil x 








Thus x4 obeys 
[18] Lixs = wxa (u = moc/h) . 


It is clear that x1, x2, and x3, which will be renamed S,, —S,, and Sy; re- 
spectively, are the components of a 3-vector. From [17], by simple factoring, 














| 5} 0 ips ipy||0 ips —ipy|| Sz 
(por+m,?c?) | —Sy = —tpz 0 1p: ips 0 —tp: 1 —S, 
Sz —tpy tp, O || tpy pz 0}; Sz 
| oO ip, ipy|| (¥ XS), 
=h|-ip. 0 ip.|| -(vxS),|, 





|—tpy —tp2 O|| (VXS)z 


so that finally 


(por +mic?)S = —hV X(¥ XS), 
or 


[19] ‘ IS iy xX (Vv XS) = —-2’S. 
c oat 
In the light of [18] and [19] it is clear that the four components of x are the 
potentials (scalar and vector) associated with a field arising from a particle of 
rest mass mp. The first-order equations have been lost sight of but they must 
be some sort of field equations in an unfamiliar form. 
Returning again to equations [15] we next eliminate x, obtaining 


0 AA* AB" 
th ae = (et Pare ) vit ve 
[20] ™ - 
* * 
ae = _ BA" 
ih i (Ent mg v2 a, V1. 
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It is to be observed that these equations are of the first order in the time 
derivative and second order in the space derivatives reminding one of Schroe- 
dinger’s equation. They may be written as 


|. ,4A AB | 
0 | fet ‘ “my 
{21] th FY, = 7 . ah? 
BA_ _, _BB 
| te 0C~*~*~S:*~<i«R 


where only the operator part has been recorded. This certainly looks like the 
Hamiltonian of a particle with both positive and negative energy states: in 
fact, except for the second-order derivatives, it reminds one very much of 
Dirac’s equation. (It is necessary to remember that each element in [21] is a 
3 by 3 matrix.) The operator in [21] is not Hermitian but the original 10 by 10 
operator is. 

The operator on the right of [21] when squared gives Ey?+c*p?, just as the 
Dirac Hamiltonian for a free particle does. The 10 by 10 matrix does not have 
this property at all. The fact that [21] behaves like the Dirac Hamiltonian 
strengthens the conviction that this equation describes the genuine particle 
aspect of the reality which [12] represents (if any such reality exists). Let us 
treat [20] as the equations for a particle, i.e assume that both y; and 2 contain 
the factor exp(—zEt/h). If E is positive and the momentum small, yo K yy, and 
the first approximation for the large component (with time dependence 
factored out) is obtained by solving 
(22 E= Ey +AA*/mo = Eo+p* /2mo, 


where only the operator has been written down. It is hardly necessary to point 
out that beyond all doubt [21] describes a particle. 

Since we have a particle to deal with, let us imagine it to be charged and 
in a magnetic field. The right thing to do seems to be to write A as a- p and 
A* as a*- p (refer to [13] to see the meaning of a) and then to subtract (e/c)A 
from each p, where A is the vector potential of the magnetic field. Then [22 
becomes 


mee natao-ta)o-fo-t0)] 


For simplicity choose A, = —}H,, Ay = }Hz, which provides a steady field 
H in the z direction. The reduction of [23] is easily carried out by using the 
following relations between the a’s: 


a,a,* = a,a,* 


ps = . S 
= 88, = 2) 

= 7  _.. 
a,a,*+a,a,* = a,a,*+a,a,* = a,ya,*+a,0,* = 0. 


(A striking similarity with the Pauli matrices for a particle of spin 3 is dis- 
played.) The final result is 
2 2ry72 


= ft . f. OH ta yt) 10." 
[24] E = Eot+ ee HM.+97 al +y°)—poll (2ia,a, ), 
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where po = eh/2moc = one Bohr magneton, and 


im. + 
(25] Qiaaz*=|i O O 
0 0 O 


This matrix has the eigenvalues 1, 0, and —1, so that the particle behaves as 
if it had unit spin and a magnetic moment of one Bohr magneton. 

Without pursuing our study any further it seems clear that [12] is an operator 
which unites in a very suggestive manner the particle and field aspects of 
what seems to be an entity of unit spin. There are numerous representations 
other than [12], and one additional representation which completely suppresses 
the particle aspect will now be presented. We subject [12] to a transformation 
by means of the matrix S, where 


ve” 3 
[26] S=S'=,|0 v2 0}. 
eae 


It is necessary to remember that S must be compatible with [12] so that, for 
example, each | is a 3 by 3 diagonal matrix, the »/2 is a 4 by 4 diagonal matrix, 
etc. The transformed matrix is 


0 M Po 
[27] imyc = | M* 0 Ne 
Po N 0 
where 
000 ->p. —py —p: 0 O 
[28] M=|0 0 0 —}py!], N= Dez 0 --—p, O|, 
000 —», 0 Pez py 0 


and fo is a diagonal 3 by 3 matrix. The 10 elements of the column matrix upon 
which [27] operates will be called, in order starting at the top, as follows: 


E,/tp, E,/in, E,/ip, —H,/p, H,/n, —H,/y, 1p, Az, Ay, A,. 


This leads to 10 simultaneous, first-order differential equations, which combine 
very naturally to give the four which follow: 


E = ab < Bite. 
H = Vv XA, 
[29] V-E= —n’¢, 


Vv XH = E/c—p?A. 
These are the equations of the field generated by a particle of spin unity and 
rest Mass Mo. 
MESON FIELD (ALTERNATIVE TREATMENT) 


It is interesting to recall that all the results achieved above stem from the 
properties of the 4-vector a which obeys (a. p)? = P*(a.p). A_ similar 
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equation holds in three dimensions for the components of the spin vector of a 
particle of unit spin. Let ¢ = ir, +jr, + kr, be a vector with components 
given by [6] and let T be any ordinary vector. Then 


[30] (e- T)* = T*(e- T). 


The expression t. T means | || T | cos @ where @ is the ‘‘angle’’ between + 
and T. Therefore one interprets [30] to mean 


[31] | ¢ |? cos*@ = cos 8. 


But | « |? = 2 so that cos@é = 0 or +1/+/2. ¢ then makes one or other of these 
angles with any vector at all. Apparently the vector a has somewhat analogous 
properties in four dimensions. In three dimensions it is possible to have a 
particle of unit spin which is really two particles of spin } bound together in 
some manner. If ¢ is the vector associated with the total spin of this system, 
then 

[32] (e- T)4 = T(<- T)?, 


as can easily be verified. Each component of ¢ has the eigenvalues 1, 0, 0, —1 
(the singlet state introduces the second zero). ¢ in this case is a 4 by 4 matrix. 
One is tempted to wonder if it would be worth looking for a 4-vector @ with 
the property that 


[33] (a- p)* = P*(a- p)’. 


The algebra of these a’s is probably very complicated and to approach the 
problem by a study of this algebra would be an arduous task. However, a 
16 by 16 matrix which obeys [33] will be written down. It has been obtained 
chiefly by guessing so the method of derivation will not be presented. That it 
has the desired property is easily verified. It is 

7 0 85 po 


‘ 0 0 0 n 

[34] imc = | 5 0 0 0 

po 1 0 0 

Each element of [34] is a 4 by 4 matrix and 

10 0 be Pez 0 -ip, O 0 
[35] 5=/0 0 0 Of}, a=|ipy 0 —pe be 
lm» @ @ & 0 -p O O 
pb, O -ip, O 0 pe O 0 


If M is the matrix in [34], then M4 = P?M?, and this is the simplest equation 
which M obeys. 

It is not easy to see, with a casual glance, that [34] leads to field equations 
but this would become obvious if one were to write down the 16 equations to 
which it gives rise. However it can be subjected to a transformation which 
puts it in a most interesting form. (The transformation matrix will not be 
given, as it is easier to check that the result satisfies M4 = P?M? than it is to 
perform the transformation.) The result is 
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le 
U : 0 
[36] Woe = Pose a bee Sey 
O° °=) ae 


where U is exactly the 10 by 10 matrix given in equation [27], the small zero 
is a single zero in the 11th row and 11th column, and 


0 Po bz Py pz 
bo O 0 0 0 
[37] V=-lp, 0 0 O O}. 
“a a 
>, ©» 6 6 6 





In the representation given by [36] the 16 functions of the column matrix 
on which [36] acts divide into three mutually exclusive groups. One function is 
zero; there is a group of 10 functions which lead once again to the set of 
equations [29]; and a group of five functions on which the operator V acts. 
It seems logical to assume that the 10 functions describe the field arising from 
particles in the triplet state and the five functions the field from the particles 
in the singlet state with zero spin. 

If the five functions on which V operates are called (in order)— x, To, 71, 
T>2, and 73, and if pz, py, pz are written as pi, p2, ps and ict, x, y, 2, aS Xo, X41, 
Xo, X3, then 





3 
oT, _ 
and 
[39] a ee ee 
© Ox, , ’ ’ ’ ® 
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NOTE 


THE EMISSION SPECTRUM OF AIF IN THE REGION 1700-2000 A 


By S. M. NAupbé anp T. J. HuGo 


In previous papers (1) the emission spectrum of AIF in the region 2000- 
8600 A has been described. Using the same excitation conditions, the emission 
spectrum in the region 1700-2000 A was photographed in the first order of a 
two-meter vacuum spectrograph with a reciprocal dispersion of 8.6 A/mm. 

The B'Z-X'Z and C'Z-—X'Z systems were both emitted strongly. The bands 
in both systems are degraded to shorter wavelengths. In the more intense 
B'Y-X'Z system, 40 bands were observed whereas in the C'Z—X'Z system 
only 10 were observed. These band-heads with their vibrational assignments 
are listed in Tables I and II. 

















TABLE I 
BAND-HEADS OF THE B!2-X'!Z sySTEM 
VA VA Va VA 
v'v”’ (cm.~) v'v”’ (cm.~) v'v"’ (cm.~) v’y"’ (cm.~) 
0,3 51941.3 1,3 52790.4 2,3 5§3624.3 4,4 54481.2 
1,4 §2027.5 2,4 52862.5 3,4 53681.7 1,0 55128.7 
2,5 §2110.2 3,5 52928.2 4,5 53731.7 2,1 55171.1 
3,6 52188.8 4,6 52988.9 5,6 53775.6 3,2 55210.5 
4,7 §2254.5 5,7 53040.7 6,7 53808.2 4,3 55240.0 
5,8 52316.1 6,8 53083.4 7,8 53828.6 5,4 55266.1 
6,9 52368.1 7,9 53113.4 0,0 54282.7 6,5 55280.3 
7,10 52407.9 8,10 §3131.9 I 54341.2 2,0 55956.6 
8,11 2431.8 0,1 53492.8 22 4393.0 3,1 55984.5 
0,2 52711.7 1,2 53561.1 3,3 54441.7 4,2 56021.3 
TABLE II 
BAND-HEADS OF THE C!2—X'!Z SYSTEM 
Vr Vh 

v'v"’ (cm.~?) v'v"’ (cm.~?) 

1,4 55571.9 0,1 56960.2 

2,5 55732.4 1,2 57104.7 

0,2 56179.4 2,3 57247.2 

1,3 56331.7 0,0 57749.8 

2,4 56482.8 tI 57884.6 


The band-heads can be represented by the following equations: 
B'>-X'Z bands: 
v, = 54249.6+[864.8(v’ +4) —7.23(v’ +4)?+0.009(v’ +4)?] 
—[799.2(v'’ +3) —4.74(0" +4)?+0.015(0" +3)] . 
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C'> —X'5 bands: 
v, = 57682.2+[934.3(v’ +3) —5.18(v’ +3)?] 
—[799.0(v"’ +3) —4.74(v"’+34)?] . 


1. Naupf£, S. M. anp Huco, T. J. Can. J. Phys. 31: 1106. 1953; 
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LETTERS TO THE EDITOR 





Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of tssue. No proof will be sent to the authors. 


On the Interpretation of Nuclear Resonance Fluorescence Experiments 


In recent years a number of observations have been made of resonance fluorescence in 
nuclei (2, 6, 7, 8, 9, 10, 11, 12). The results have been interpreted in terms of the basic resonance 
formula relating an absorption cross section ¢ to a transition probability (times) I: 


n2 r 
i) ° = £8, (E—Ey* +072) 
It has been assumed that the data so interpreted yield directly the partial mean life of the 
transition for gamma radiation (r = #/T) even in the presence of appreciable branching and/or 


internal conversion. ; 
It is the purpose of the present note to point out that this interpretation is not valid. In the 
presence of branching or internal conversion, the absorption cross section for the resonant 


absorption and subsequent re-emission of a gamma ray is: 
2 we: dle iacsiniaa lea 
[2] o = 885 (E—Eo)* +y/2)8 


where y: is the gamma ray transition probability for the absorption (and re-emission) transi- 


tion, while y = > vit): a is the total width of the upper level including all possible modes 
of de-excitation to all lower levels. y; is the gamma transition probability for the 7th gamma 
ray originating from the upper level, while a; is the corresponding internal conversion transition 
probability. 

In the resonance fluorescence experiments it is the integral of the cross section over energy 

that matters. Therefore from Equation [1] we find §«¢dE~ TI, while from Equation [2] we 
obtain §{o dE~ 7:°/y. If we denote the ‘‘gamma ray partial mean life” of the upper level as 
obtained in the past from Equation [1] as (ry):, it follows that the correct gamma ray partial 
mean life is given by 
[3] Ty = (y1/y) (ty) 1- 
In the simple case where the upper level is the first excited state then (y:/y) = 1/(1+etotai) 
where aiotai is the total internal conversion coefficient. If the mean life of the level (including 
conversion) is desired, ry given by Equation [3] must be multiplied again by (71/7). If the upper 
level is not the first excited state the ratio (y:/y) will involve branching ratios as well as 
conversion coefficients but the analysis will be the same. 

The experimental data available at present are not very precise, but the modification implied 
by Equations [2] and [3] tends to improve the agreement between the ry values obtained from 
resonance fluorescence and those obtained from delayed coincidence experiments (1, 3, 4). 
An example is the 209 kev. transition in Hg!**. Metzger (8) quotes a ‘‘mean life’”’ based on 
Equation [1] of (ry): = (83+1) X10~-” sec. The de-excitation of the 209 kev. level involves 
both branching and conversion (13). The experimental value for (y:i/y) is approximately 0.5. 
Thus the correct gamma ray partial mean life is 7, = (1.5+0.5) X10~™ sec., and the mean 
life of the 209 kev. level is 7 = (7.5+2.5) X10-" sec. This can be compared with Graham's 
delayed coincidence measurement (5), 7 = (10+6) X10-" sec. It is evident that the large 
experimental uncertainties in this and other cases preclude accurate comparisons. As the data 
become more accurate, however, the difference between the “‘lifetime’’ obtained from Equation 
[1] and the correct 7, obtained from Equation [2] will be significant for those transitions in 
which branching and/or conversion are appreciable. 


The author wishes to thank Dr. R. E. Bell for drawing his attention to this problem and for 
helpful discussion about it. 
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